EXACT SOLUTION OF THE SIX- VERTEX MODEL WITH DOMAIN 
WALL BOUNDARY CONDITIONS. ANTIFERROELECTRIC PHASE 



PAVEL BLEHER AND KARL LIECHTY 

Abstract. We obtain the large n asymptotics of the partition function of the six- vertex 
model with domain wall boundary conditions in the antiferroelectric phase region, with the 
weights a — sinh(7 — t), b ~ sinh(7 + t), c~ sinh(27), \t\ < 7. We prove the conjecture of 
Zinn- Justin, that as n ^ 00, Z„ ~ C'di{nuo)F'^ [1 + 0(n^^)], where w and F are given by 
explicit expressions in 7 and i, and "dAiz) is the Jacobi theta function. The proof is based 
on the Riemann-Hilbert approach to the large n asymptotic expansion of the underlying 
discrete orthogonal polynomials and on the Dcift-Zhou nonlinear steepest descent method. 



1. Introduction and formulation of the main result 

1.1. Definition of the model. The six-vertex model, or the model of two-dimensional ice, 
is stated on a square n x n lattice with arrows on edges. The arrows obey the rule that at 
every vertex there are two arrows pointing in and two arrows pointing out. Such a rule is 
sometimes called the ice-rule. There are only six possible configurations of arrows at each 
vertex, hence the name of the model, see Fig. 1. 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



Figure 1. The six arrow configurations allowed at a vertex. 
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We will consider the domain wall boundary conditions (DWBC), in which the arrows on 
the upper and lower boundaries point into the square, and the ones on the left and right 
boundaries point out. One possible configuration with DWBC on the 4x4 lattice is shown 
on Fig. 2. 



Figure 2. An example of 4 x 4 configuration with DWBC. 

For each possible vertex state we assign a weight Wi, i — 1, . . . , 6, and define, as usual, 
the partition function as a sum over all possible arrow configurations of the product of the 
vertex weights, 

Zn= E ^(^)' ^(^)= n^*(-;<^)=n^f'^''^ (1-1) 

arrow configurations a xEVn i=l 

where Vn is the nxn set of vertices, t{x; a) G {1, . . . , 6} is the vertex-type of configuration a at 
vertex x according to Fig. 1, and Ni{a) is the number of vertices of type i in the configuration 
a. The sum is taken over all possible configurations obeying the given boundary condition. 
The Gibbs measure is defined then as 

/-„W = ^. (1.2) 

Our main goal is to obtain the large n asymptotics of the partition function Z^- 

The six- vertex model has six parameters: the weights Wj. By using some conservation 
laws it can be reduced to only two parameters (see, e.g., [14], [7], and [1]). Namely, we have 
that 

Z„(wi, W2, W3, W4, W5, We) = C{n)Zn{a, a, b, b, c, c), (1.3) 

and 

Ai„(cr; wi, W2, W3, W4, W5, Wq) = Ai„(cr; a, a, b, b, c, c), (1.4) 

where 

a = ^/vhW^, b = y/W3W4, C = y/w^WQ, (1.5) 
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and 



Furthermore, 



and 



We, 



Zn{a, a, b, b, c, c) = c" Z„ 



a a b b 

) ) ~) ~5 1) 1 

c c c c 



(1.6) 
(1.7) 



/in(o-; a, a, ^ b, c, c) = /x„ ( a; -, -, -, -, 1, 1 ) , 

c c c c 



a b 

c 



SO that a general weight reduces to the two parameters, ^, 

1.2. Exact solution of the six- vertex model for a finite n. Introduce the parameter 

+ 6^ - 



2ab 



(1.9) 



The phase diagram of the six-vertex model consists of the three phase regions: the ferro- 
electric phase region, A > 1; the antifcrroelectric phase region, A < — 1; and the disordered 
phase region, — 1 < A < 1. Observe that |a — 6| > c in the ferroelectric phase region and 
c > a -|- 6 in the antifcrroelectric phase region, while in the disordered phase region a, b, c 
satisfy the triangle inequalities. 








1 



a/c 



Figure 3. The phase diagram of the model, where F, AF and D mark fer- 
roelectric, antifcrroelectric, and disordered phases, respectively. The circular 
arc corresponds to the so-called free fermion hne, when A = 0, and the three 
dots correspond to 1-, 2-, and 3-enumeration of alternating sign matrices. 



In the three phase regions we parameterize the weights in the standard way: for the ferro- 
electric phase. 



a = sinh(t — 7), 6 = sinh(t -|- 7), c = sinh(2|7|), 0<|7|<t; 



;i.io) 
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for the antiferroelectric phase, 

a — smh.{'y — t) , b — smh.{'y + t) , c = smh(27), |t|<7; (l-H) 

and for the disordered phase 

a = sin(7 — i), 6 = sin(7 + t), c = sin(27), |t| < 7. (1-12) 

The phase diagram of the six- vertex model is shown on Fig. 3. The phase diagram and 
the Bethe Ansatz solution of the six-vertex model for periodic and antiperiodic boundary 
conditions are thoroughly discussed in the works of Lieb [21]-[24], Lieb, Wu [25], Sutherland 
[28], Baxter [2], Batchelor, Baxter, O'Rourke, Yung [3]. See also the work of Wu, Lin [29], 
in which the Pfaffian solution for the six- vertex model with periodic boundary conditions is 
obtained on the free fermion line, A = 0. Brascamp, Kunz and Wu [6] prove the equality of 
the free energy with periodic and free boundary conditions under various conditions on a, b, c, 
and also they prove the existence of the spontaneous staggered polarization for sufficiently 
small values of the parameters a/c and b/c. 

As concerns the six- vertex model with DWBC, it was noticed by Kuperberg [20], that on 
the diagonal, 

^ = ^ = (1.13) 
c c 

the six-vertex model with DWBC is equivalent to the s-enumeration of alternating sign 
matrices (ASM), in which the weight of each such matrix is equal to s^^, where is the 
number of (— l)'s in the matrix and s = The exact solution for a finite n is known for 1-, 
2-, and 3-enumerations of ASMs, sec the works by Kuperberg [20] and Colomo-Pronko [11] 
for a solution based on the Izergin-Korepin formula. A fascinating story of the discovery of 
the ASM formula is presented in the book [10] of Bressoud. On the free fermion line, 7 = f , 
the partition function of the six-vertex model with DWBC has a very simple form: Z„ — 1. 
For a nice short proof of this formula see the work [11] of Colomo-Pronko. 

In this paper we will discuss the antiferroelectric phase region, and we will use parameter- 
ization (1.11). The parameter A in the antiferroelectric phase region reduces to 

A = -cosh(27). (1.14) 

The six- vertex model with DWBC was introduced by Korepin in [17], who derived impor- 
tant recursion relations for the partition function of the model. These recursion relations 
have been solved by Izergin [15], and this lead to a beautiful determinantal formula for the 
partition function with DWBC. A detailed proof of this formula, usually called the Izergin- 
Korepin formula, and its generalizations are given in the paper of Izergin, Coker, and Korepin 
[16]. When the weights are parameterized according to (1-11), the Izergin-Korepin formula 
is 

^ ^ [sinh(7 - t) sinh(7 + ^)]"V^ 
where t„ is the Hankel determinant. 
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and 

m = . ""^'^1, — (1.17) 

^ smh(7 - 1) smh(7 + 1) ^ ^ 

An elegant derivation of the Izergin-Korepin formula from the Yang-Baxter equations is 
given in the papers of Korepin, Zinn- Justin [18] and Kuperberg [20] (see also the book of 
Bressoud [10]). 

One of the applications of the determinantal formula is that it implies that the function 
Tn solves the Toda equation 

d 

T"n< - = ^n+lTn-l, U > 1, = ^ ' (1-18) 

cf. the work of Sogo, [27]. The Toda equation was used by Korepin and Zinn-Justin [18] 
to derive the free energy of the six-vertex model with DWBC, assuming some ansatz on the 
behavior of subdominant terms in the large n asymptotics of the free energy. 

Another application of the Izergin-Korepin formula is that t„ can be expressed in terms 
of a partition function of a random matrix model and also in terms of related orthogonal 
polynomials, see the paper [31] of Zinn-Justin. In the antiferroelectric phase the expression 
in terms of orthogonal polynomials can be obtained as follows. For the evaluation of the 
Hankel determinant, let us write (j){t) in the form of the Laplace transform of a discrete 
measure, 

0(t) = . ^^^y^^^, , = 2 V e--^l^l. (1.19) 

^ ^ sinh(7 - 1) sinh(7 + t) ^ ^ ^ 



l=—oo 



Then 



where 



Tn-— E A(/)^n^'''"'*' (1-20) 



nl 

ll,...,ln =—0O i=l 



A(/)= n (^^•-^^) (1-21) 

l<i<j<n 

is the Vandermonde determinant, see Appendix A in the end of the paper. 

Introduce now discrete monic polynomials Pj{x) — + . . . orthogonal on the set Z with 
respect to the weight, 

«;(Z) = e2«-2Tl'l, (1.22) 

so that 



J2 Pj{l)Pk{l)w{l) = hSjk. (1.23) 

l=—oo 

Then it follows from (1.20) that 

n-l 

Tn^2^"l[h, (1.24) 



k=0 

see Appendix B in the end of the paper. 
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1.3. Rescaling of the weight. Set 

A„ = ^, x = /A„, wn{x) = e-<\^\-^^\ C = -, (1-25) 
n 7 

and 

P„,(x) = A^P, (^^) . (1.26) 



L^^^Jx^^, A;eZ|>. (1.27) 



Consider also the lattice 

Ln — — - 

I ^ 

Then from (1.23) we obtain that the monic polynomials P„fc(x) satisfy the orthogonality 
condition, 

^ Pnj{x)Pnk{x)Wn{x)/^n = Kk^jk, Kk = hkA.f+\ (1.28) 

We can then combine equations (1.15), (1.24), and (1.28) to obtain 

^-=(^) Po(^' " = b = smh{^ + t). (1.29) 

For what follows we will need to extend the weight Wn{x) to the complex plane. We do so 
by defining Wn{z) on the complex plane as 

w^{z) = e-"^(^) (1.30) 

where 

{z — (^z when Rez>0, 

-z-Cz when Rez<0, ^^"^^^ 
so that V{z), and thus Wn{z), is two- valued on the imaginary axis. 

1.4. Main result: Asymptotics of the partition function. This work is a continuation 
of the work [4] of the first author with Vladimir Fokin and [7], [8] of the authors of the 
present work. In [4] the authors obtain the large n asymptotics of the partition function 
in the disordered phase. They prove the conjecture of Paul Zinn- Justin [31] that the large n 
asymptotics of Z„ in the disordered phase has the following form: For some £ > 0, 

Zn = Cn'^F"'[l + 0{n-% (1.32) 
Furthermore, they find the exact value of the exponent k, 

- ^ 1.33) 
12 37r(7r-27) ^ ^ 

The value of F in the disordered phase is given by the formula, 

F = -^^, a = sin(7-i), 6 = sin(7 + i) , (1.34) 
27 cos 2^ 

in parameterization (1.12). 
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In the work [7] we obtain the following large n asymptotic formula for Z„ in the ferroelectric 
phase region: For any £ > 0, 



l + O 



(e-»") 



(1.35) 



where C = 1 — e~^^, G = e'''"*, and F = sinh(t + 7) in parameterization (1.10). 

Finally, in the work [8] we obtain the following large n asymptotic formula for Z„ on the 
border line between the ferroelectric and disordered phase regions: 

(a, a, a + 1, a + 1, 1, 1) = Cn'^G^F^'fl + 0{n-^/^)] , 

where C > 0, 



K — 



G — exp 



-c 



a + 1, 



;i.36) 
:i.37) 



and C is the Riemann zeta-function. 

In the present paper we obtain the large n asymptotic formula for Z„ in the antiferroelectric 
phase region. The formulation of the main result of the present paper and the proofs involve 
the Jacobi theta functions. Let us review their definition and basic properties. 

There are four Jacobi theta functions: 



2 ^(-l)'^g("+5)' sin ((2n + 1)^) 

n=0 

00 

2^g("+^)'cos((2n + l)z) 



n=0 



(1.38) 



^i(^) 

^^{z) = 1 + 2 ^ g"' cos(2n^) 

n=l 

00 

= 1 + 2 J](-l) cos(2n^) 

n=l 

where q is the elliptic nome. We will assume that 1 > g > 0. Fig. 4 shows the graphs of 
t?!, ■(?2 (left figure) and -i^s, ■(?4 (right figure) on the interval [0,7r] for q = 0.5. Observe that 
1?!, '^4 are increasing on [0, ^] while '1)2, "^s are decreasing on this interval. 

The Jacobi theta functions satisfy the following periodicity conditions: 

t9i(z + vr) = -Mz), + ^r) = -e-''''q-^Mz): 



+ tt) = -Mz), + vrr) = e-^''q-^M^), 

+ tt) = M^), + ttt) = e-2-g-i^3(z), 

+ tt) = M^), + ttt) = -e-^^'q-^Mz), 

where r is a pure imaginary number related to q by the equation, 

q = e'^'r 

The theta functions also satisfy the symmetry conditions: 

M-z) = -Mz), M-z) = Mz), M-z) = Mz), M-z) = M^), 



;i.39) 



(1.40) 
(1.41) 
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and the equations, 

^i(^) = ^2 - I) , ^3(^) =^4(^ + 1), ^i(z) = -ie-+^^4(^+ (1-42) 

The only zeroes of the theta functions are 

^i(0) = 0, ^2(f)=0, ^3(| + y)=0, ^4(y)=0, (1.43) 

and their shifts by rmr + mrr; m,n & Z. There are many non-trivial identities satisfied by 
the theta functions. A list of those identities used in this paper is given in Appendix E. 

In the antiferroelectric phase region we use parameterization (1.11) with two parameters 
t, 7 such that \t\ < 7. In what follows we will also use the following parameters: 

C=-e(-l,l), a; = ^^il±^e(0,7r). (1.44) 
7 2 

The elliptic nome for all Jacobi theta functions in this paper will be equal to 

g = e 27. (1-45) 
Our main result in the present paper is the following asymptotic formula for Z„: 
Theorem 1.1. As n ^ 00, 

Zn = C^4{nuj)F''\l + 0{n-^)), (1.46) 
where C > is a constant, and 

7rsinh(7-t) sinh(7 + t)^;(0) 

^ ■ ^'-''^ 

The asymptotic formula (1.46) proves the conjecture of Zinn- Justin in [31]. The proof of 
Theorem 1.1 will be based on the Riemann-Hilbert approach to discrete orthogonal polyno- 
mials. An important first step in this approach is a construction of the equilibrium measure. 



exact solution of the six- vertex model 
2. Equilibrium measure 



2.1. Heuristic motivation and definition of the equilibrium measure. If we scale 
the variables in (1.20) as fii — then we can rewrite formula (1.20) as 



2"' 



= E e--'^-M^ (2.1) 

where 



T). 



1 

du^{x)^-J2S{x-fXj), (2.2) 

and 

H{u) =11 log ^ du{x)du{y) + / {\x\ -Cx)du{x). (2.3) 

where all integrals are over M. 

Due to the factor (— n^) in the exponent of (2.1), we expect the sum, in the large n limit, 
to be focused in a neighborhood of a global minimum of the functional H. Clearly, we have 
that is a probability measure and 

i^/i(a, b) < ^ ^ for any — oo<a<6<oo, (2.4) 

because in (2.2), fXj e and fXj ^ ji^ '^^ j ^ k. With these constraints in mind, we define 

Eo = \niH{v) (2.5) 

V 

where the infimum is taken over all probability measures satisfying (2.4). It is possible to 
prove that there exists a unique minimizer i/q; so that 

Eo = Hiuo), (2.6) 

see, e.g., the works of Saff and Totik [26], Dragnev and Saff [13] and Kuijlaars [19]. Fur- 
thermore, uq has support on a finite number of intervals, and is absolutely continuous with 
respect to the Lebesgue measure. The minimizer i/q is called the equilibrium measure. 

Denote the density function of the equilibrium measure as p{x), and its resolvent as a;, so 
we have 

— = p(x), -iz) = J^, (2.7) 

and 

p(x) = —(u;(x -iO) -u;(x + iO)). (2.8) 

The structure of the equilibrium measure vq is studied in the paper of Zinn- Justin [31] , who 
shows that has support on an interval [a, /3], with a saturated region [a', /3'] in which 

pix) = ^, xe[a',(3'], (2.9) 

and two unsaturated regions, [a, a'] and [/?',/?], in which 

0<p{x)<^, X e {a,a')Li{P',P), (2.10) 
27 
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see Fig. 5. We also have that 



a < a' <0 < (3' < (3, 



so that the origin, which is a singular point of the potential V{x) 
saturated region [a', P'] . 



(2.11) 

\x\ — Cx, lies inside the 





P(x) 


• ^ 




a a' 


P' p 



Figure 5. The equilibrium density function p{x). 

The measure z/q is uniquely determined by the Euler-Lagrange variational conditions 

{= / for a; G [a, a'] U [/?',/?], 
>/ for ,TG[a',/3'], (2.12) 
<l for x^{a,p\, 

where I is the Lagrange multiplier. The Euler-Lagrange variational conditions imply 

a;(x-iO) +a;(x + iO) = -C + sgn(x) for x e [a, a'] U [/?', /?], (2.13) 
whereas in the saturated region, we have 

(2.14) 



p(x) = — (a;(x-iO) -a;(x + iO)) = — for x e [a', 
27ri 27 



Now we will give a detailed description of the equilibrium measure. We begin with explicit 
formulae for the end-points of the support of the equilibrium measure. 

2.2. Explicit formulae for the end-points. They are given in the following proposition. 
Proposition 2.1. The end-points of the support of the equilibrium measure uq are equal to 

(2.15) 



a — — TT ■ 



= -TT 



mi 



/3= -TT 



mi 
m) 
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The differences between the end-points are equal to 



11 



and 



3 , ,,2.o /3(f)^4(f) 



(2.16) 



i ^l(f)^3(f) 



(2.17) 



^l(!)^3(f)' 

Finally, we have the Zinn- Justin formula for the centroid of the end-points, 

a + a' + 13' + (3 _ vr^^'slf) 



4 2^2(f)' 
For a proof of Proposition 2.1 see the next section. 



(2.18) 



2.3. Equilibrium density function. The equihbrium density function is described in the 
following proposition. 

Proposition 2.2. The equilibrium density function p{x) is given by the formulae, 

( 1 r dx' 



Ja ^/{x' - a){a' - x'){(5' - x'){(5 - x') 



, a < X < a' , 



p{x) = < 



— , a'<x<P', 
27 



1 



(2.19) 



d.r' 



Also, 



[ T^Jx ^/{x' - a){x' - a'){x' - (5'){(5 - x') 

+ c 



, 13' <x<l3. 



I p{x) dx — — 
Jo 



The resolvent u}{z) of the equilibrium measure is given as 

I dz' 

uj{z) 



^{z' - a) {z' - a') {z' - (3') {z' - p) ' 



(2.20) 



(2.21) 



where integration takes place on the sheet of \/ R{z') = ^/ [z! — a) {z' — a') {z' — (3') {z' — j3) 
for which ■\jR{z') > for z' > (3, with cuts on [a, a'] and [13' , 



For a proof of this proposition see the next section. 
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2.4. g'-function. Define the g'-function on C \ [—00, P] as 



g{z)^ I \og{z - x)dvo{x) (2.22) 

J a 



where we take the principal branch for logarithm. 
Properties of g{z): 

(1) g{z) is analytic in C \ (— oo,/3]. 

(2) For large z, 



g{z) = \ogz-Y,-,. gj= -duoix). (2.23) 

(3) 

g'iz)=u;{z). (2.24) 
(4) From the first relation in (2.12) wc have that 

g+{x) + g-{x) = \x\ - (x + I for x e [a,a']U [j3' , jS], (2.25) 

where (?+ and refer to the limiting values of g from the upper and lower half-planes, 
respectively. By differentiating this equation we obtain that 

uj+{x)+uj-{x)^g'_^.{x)+g'_{x)^sgnx-C for x e [a,a'][J [p' , (3], (2.26) 
Consider the function 

fix) = g+{x) + g_{x) - (Ixl - + 0- (2-27) 
We have from (2.25), (2.26) that 

f{x) = f{x) = for x = a, a', (3', (3, (2.28) 
and from (2.21) that 

/» = --= ] for xe{-oo,a)U{a',f3')U{f3,oo). (2.29) 
^/ [x — a)[x — a'){x — p ){x — p) 

Since 

f"{x) < for xe (-00, a) U (/3, 00). (2.30) 

and 

f"{x) > for X e {a', (3'), x ^ 0, (2.31) 

we obtain that 





r = x — + / 


for 


X e [a, a'] U 




g+{x)^g^{x) 1 


>\x\-C,x + l 


for 


X e (q;',/3'), 


(2.32) 




[ < x - Cx + z 


for 


X e R \ [a,/?]. 





(5) Equation (2.22) imphes that the function 

G{x) = g+{x)-g_{x) (2.33) 
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is pure imaginary for all real x, and 

' 27ri for — oo < X < a, 

2m — 2TTi / p{s) ds for a < x < a, 

J a 

'1 + C x' 



G{x) 



2m 



27 



for a' <x< f3', 



(2.34) 



2m / p{s) ds for P' <x<P, 

J X 



for f3 < X < oo. 
From (2.32) and (2.34) we obtain that 

\x\ — (x + I ± 



2g^{x) = < 
(6) Also, from (2.34) 



rx 

2m — 2m I p[s) ds 

J a 

\x\-Cx + l± 2m / p{s) ds for P' < x < p. 

^ J X 



for a < X < a', 



(2.35) 



dG{x + iy) 



dy 



2np{x) > 0, a; e {a,P). 



(2.36) 



Observe that from (2.25) we have that 

G{x)^2g^{x)-V{x)-l^-[2g_{x)-V{x)-l], x e [a, a'] U (2.37) 
where V{x) = \x\ — (x. 

2.5. Evaluation of the Lagrange multiplier /. Wc have the following proposition. 
Proposition 2.3. The Lagrange multiplier I solves the equation, 



i 

62 



(2.38) 



For a proof of this proposition see the next section. 

3. Proof of Propositions 2.1, 2.2, and 2.3 

Proof of Proposition 2.1. Following Zinn- Justin [31], we make the following elliptic 
change of variables: 

^ r , (3.1) 

^ \/ [z' - a)[z' - a')[z' - P')[z' - P) 



where integration takes place on the sheet on ^y R(z') specified in Proposition 2.2. To 
understand this integral in terms of the usual elliptic integrals, we first make the change of 
variables 

, (P - z'W - a) 



14 PAVEL BLEHER AND KARL LIECHTY 

SO that 

^ {(5 - a)v - {(5' - a) ' ^ ' ' 

Note that = 0, v{l3') — oo, and v{a) — 1. When we substitute v into equation (3.1), 
we have 

= ^ / I ' (3-4) 

where 



•v{z) 




dv 


\ 


Iv{v - 




'(P- 


a){P' 


-a') 


(/?'- 


-a){(3 


-a') ' 



(3-5) 

We next take v' = ^/v, obtaining 

u{z)= / , (3.6) 



which corresponds to ^/v{z) = sn{u,k), so that 

if3 - z){f3' - a) 



{f3'-z){P-a) 



sn (m), sn(M) = sn(M, /c). (3.7) 



Notice that u maps the upper z-plane conformally and bijectively onto the rectangle [0, K] x 
[0, iK'] , and the lower 2;-plane conformally and bijectively onto the rectangle [0, K] x [—iK', 0] , 
where 

K — uia) — I — = and 

Jo Jil-v'^){l-kW 

^\ ' (3.8) 

/" *: dv' 

K' = -iu{(3') = / _ 

are the usual complete integrals of the first kind. More specifically (see Fig. 6), 

(1) The upper (resp. lower) cusp of the interval [13' , (3] is mapped onto the interval [0, iK'\ 
(resp. [0,-iK']). 

(2) The upper (resp. lower) cusp of the interval [a, a'] is mapped onto the interval 
[K,K + iK'] (resp. [K,K-iK'\). 

(3) The interval [a', f3'] is mapped onto the inteval [iK' , K+iK'] or the inteval [—iK', K— 
iK'], depending on the path of integration. 

(4) The remaining part of the real axis, [—00,0;] U 00], is mapped onto the interval 
[0, K], with ■u(oo) = = Moo- 

We will denote the rectangle [0, K] x [—iK' , iK'] as R, the fundamental domain of the function 
z{u). We can now define 

i:j{u) = Lj{z{u)) for ue R. (3.9) 
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u(P')=iK' a 



u(P)=0 



u(P')=-iK' 



u(a')=K+iK' 




u(a)=K 



u(a')=K-iK' 



Figure 6. The i^-plane. Here — u^o = u{oo) and a — u^-\- iK' . 



The Euler-Lagrange equation (2.13) and the equation (2.14) then become 

a;(M) + Cj{-u) = 1 - C for e [-iK\ iK'] 

u{u) + u{-u + 2K)^-1-C for u e [K - iK', K + iK'] 

iTT 



Lj{u + 2iK') -uj{u) = -- 



7 



for M e [-iK' , K - iK']. 



(3.10) 



The function uj{z) is analytic in C\[a, /?], but can be analytically continued from either above 
or below through any of the cuts [a, a'], [a',/?'], and These analytic continuations in 

the 2;-plane give an analytic continuation of a; in the rt-plane into a neighborhood of R, which 
can then be continued by equations (3.10) to the entire w-plane. We therefore have that u 
is analytic and satisfies equations (3.10) throughout the u-plane. The first two equations of 
(3.10) can be combined as 



(3.11) 



It therefore follows that a; is a linear function of u, as its derivative is a doubly periodic 
entire function. We also know from the fact that uj{z) ~ - at infinity that 



^J{(5' -a){(5-a') 



[U - Uoo) + 0{U - Uoof 



(3.12) 



in some neighborhood of Mqo, where Mqo is the image of infinity under the map u{z). It thus 
follows from (3.10), (3.11), and (3.12) that 



LJ[U) = — - 



1 

K' 



[U - Ur_ 



(3.13) 
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and that 



K' _ t: 
^ ~ 27 

^/{(3' -a){(3-a')^2K (3.15) 

Uoo ^ 1 - C 
K ~ 2 ' 



(3.14) 



(3.16) 



From (3.7) we obtain that 
This imphes that 



1 - sn (Moo) = 1 - 



^-^=s^\u^). (3.17) 

p ~ a 

13' -a 13-13' 



13 — a (3 — a ^ 
dn.\uoo) = 1 - k^sn^iuoo) = 1 - 



2 2, ^ , {(3-a){(3' -a'){(3' -a) (3 - (3' 



(3.18) 



{(3' - a){(3 - a') {(3 - a) (3 -a'' 
Prom equations (3.15), (3.17), (3.18) we obtain the distances between the turning points 



in terms of u. 



oo- 



sn(Moo)cn(Moo) 



P-a'^ 2K (3.19) 

sn(Moo)an(Moo) 



sn(uoo) 

The functions sn, cn, and dn are expressed in terms of Jacobi theta functions as follows (see, 
eg. [30]), 

MO) MSI MO) MSI MO) MSI con^ 

By (3.14), the half-period ratio r and the elliptic nome q of the theta functions are 

iK' in , -nK' -TT^ 

T = = — and q = eK = e 27 . (3.21) 

If we take into account the fact that 

2K 

MOf = — (3.22) 

TT 

and equation (3.16), we can write equations for the distances between the turning points 
that involve the original parameters only: 

(3.23) 



M%)M%) ' " ' ' M%)Mi) 



'M%)M%y 
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giving (2.17). These equations determine the end-points a, a', /3 up to a shift. To fix the 
shift we use the equation (2.12) at the points a' and (3' to obtain 

/ {uj{z + iQ)+uj{z-iQ))dz = {l-Q(5' + {l + Qa'. (3.24) 

Ja' 

Writing this integral in terms of u gives 

/ —{u-u^)r\u)du+ —{u-u^)r'{u)du = {l-Q(5' + {l + Qa' (3.25) 
where 

((3 - a)sn^iu) - iP' - a) ' ^ _ r [u) - ^^r{u). {6.2b) 



Note that r{±iK') = j3' and r(JC ± iK') = a' . Integrating by parts gives 
((K ^ B'v.A - I 

K 



{{K-u^)a' + P'u^)- / -^du- / -S^du = {l-QP'+{l + Oa' (3.27) 
or equivalently 

pK+iK' nK-iK' 

/ r{u)du+ / r{u)du^Q. (3.28) 

JiK' J-iK' 

We can evaluate these integrals by first writing r{u) in the form 

r M = + . 3.29 

and using the functions 

The addition formulae for the sn and Z functions arc (sec [30]): 

, , , sn(M)cn(a)dn(a) ± sn(a)cn(M)dn(M) 

snfii ±a) — — — ' o/ N O/ N — , 

1 - Psn2(a)sn2(M) (3.31) 

Z{u ± a) = Z{u) ± Z{a) =F A;^sn(ti)sn(a)sn(ti ± a). 

Thus we have 

Z{u — a) — Z[u + a) + 2Z{a) = A:^sn(-u)sn(a) (sn(-u + a) + sn(-u — a)) 

A;^sn(?7,)sn(a) [2sn(M)cn(a)dn(a)] 
" 1 - fc2sn2(a)sn2(M) (3.32) 

2A;2sn(a)cn(a)dn(a)sn2(ti) 
1 — k'^sn'^{a)sn'^{u) 
We also have the half- and quarter-period identities 

sn(u + iK') = — , cn(u + iK') = Z!^!^ , dn(u + iK') = -z^^. (3.33) 
A;sn(M) k sr[u) sn[u) 
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In particular, notice that ^^^^ ^ = ksn{u^ + iK'). Using the addition formulae (3.32) we 
can write r{u) as 

r{u) =13+ f ^,, , , , „, - a) - Z{u + a) + 2Z{a)) , (3.34) 



2A;2sn(a)cn(a)dn(a)sn2(woo) 
where a = u^o + iK', (see Fig. 6). Prom (3.33) and (3.19), it follows that 

Thus we can write (3.34) as 

r{u) = 13-K [Z{u - Moo - iK') - Z{u + u^o + iK') + 2Z{u^ + iK')] . (3.36) 

If we write u = x + iK' in the first integral of (3.28), and u = x — iK' in the second, we 
obtain 



/•ft 

/ [2/3 - AKZ{u^ + iK') - K[Z{x ~ u^) - Z{x + + 2iK') 
Jo 

+Z{x — Uoo — 2iK') — Z{x + Moo)]] o?a; = 
Prom the periodic properties of ■j94, it follows that 



(3.37) 



Z{u±2iK') = Z{u) + —, (3.38) 

K 



so we can write (3.37) as 

^ [2{i - 4:KZ{uoo + iK') - 2m + 2K [Z{x + Moo) - Z{x - u^)] ^ rfa; = 



(3.39) 



This equation is readily integrated, as Z is the logarithmic derivative of the © function. 
Integrating gives 

K 

«{ .T -I- 7;.._ I 





{2(3 -AKZ{u^ + iK') -27ri)x + 2K log ^ 



Q{x '^ooj_ix=o (3 40) 

= 2Kf3 - 4K'Z{u^ + iK') - 2Kni + 2K\og /^ Qj^ + M ©(-m ) \ 

\e{K-u^) e(Moo) / 

The logarithmic term in this equation is zero due to the evenness and periodicity (period 
2K) of the function and the fact that the relevant term in the integration is real on the 
entire contour of integration. Thus we have that 

/3 = 2KZ{u^ + iK') + ni. (3.41) 

Prom (1.42), we can deduce that 

Ziu^ + .K')=-^m.^ (3.42) 



and write (3.41) as 
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This equation, together with equations (3.23) and (3.21), determine the end-points, a, a', f3. 
In fact, similar to (3.43) we have the following explicit formulae for the other end-points: 



a — — TT 



2 ^ 



a — — TT ■ 



(3.44) 



^l(f)' ^ "^4(f)' " "^3(f)- 

This follows from (3.23), (3.43), and the identities (E.2). Similarly, in addition to the 
formulae (3.23) for distances between turning points, we get (2.16). 

Proof of the formula of Zinn- Justin for ^±^l±SL±Ii_ vVe immediately have from (3.44) that 



a + a' + P' + P 



TT 
"I 



M%) M%) U%) M%). 



Prom the identity (E.6), we can deduce that 

i?;(2z) _ 1 

thus we have 



a + a' + /3' + /3 



Mz) Hz) M^) M^). 

^ TT ^Uf + f) ^ 7r^^,(f) 

2'di{u;) 



(3.45) 



(3.46) 



(3.47) 



2^i(f + f) 2^,(f) 

Proof of Proposition 2.2. From equations (3.13), (3.1), and (3.15) wc obtain formula 
(2.21), cf. [31]. From formula (2.21) and equations (2.8), (2.14) wc obtain that the equi- 
librium density function p{x) is given by formulae (2.19). We are left to prove formula 
(2.20). 

By (3.1), (3.15), and (2.19), on the interval 

p{x) = j^u+{x) for xe[i3',i3]. (3.48) 



zKtt 



It follows that 



J0' 



p{x)dx 



iKtt 



-3 



u+{x)dx 



ur'{u)du, 



(3.49) 



K' 



where r{u) is defined in (3.26). If we use equation (3.36), together with formula (3.41), we 
can write r{u) as 

r{u) — m — K [Z{u — — iK') — Z{u + u^o + iK')\ 
Integrating (3.49) by parts, we get 



JiK> 



ur'{u)du = -p'iK' - ttK' - K 



log 



e{-uoo - iK')eiuoo + 2iK'] 

Q{-u^)e{u^ + iK') 



Using the fact that © is an even function and the identity 



Q{u + 2iK') = e'^e-'^^e—Q{u), 



we can write (3.51) as 



ur'{u)du = -P'iK' - nK' + K [{7: + tt— - 



iK' 



K' 



K 



(3.50) 
(3.51) 
(3.52) 

(3.53) 



i(/s:7r - 13' K' 
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Remark: There is some question here as to which branch of the logarithm to take, but it is 
clear that we have chosen the correct branch, as it is the only one that gives < jf p{x)dx < 



1. 



Thus, from (3.49) and (3.53), we have 



J0' 



p{x)dx 



iKir 



i{KT: - (5'K' - T^u^ 



1 - 



ttK K 



hence by (2.19), 



['^ Z"^' /"^ (3' 1 - C 
p{x)dx = p{x)dx + j^^ p{x)dx =2^ + 1-2^ 2~ 

which proves formula (2.20). 

Proof of Proposition 2.3. By taking x = P we obtain from (2.32) that 



27 2 ' 



-C 1 + C 



We also have that 



hence 



/ = 2g{f3) - V{f3) = 2g{f3) - (1 - QP- 



lim [^(^) -log^] = 0, 

A— »oo 



I - -2 lim [giA) - g{/3) - log A] - (1 - ()/? 

A— +00 



= -2 hm 



uj(z) dz — log A 
Writing this integral in terms of u {so z — r{u)) gives 



(1-0/3. 



Z = 2 hm 



2 lim 



I — {u - u^y {u) du + log A 
Jo ^ 



-(1-0/5 



K 



{u — Uoo)r'{u) du + log r(i?) 



-(1-0/5, 



where A — r{B). Integrating by parts gives 

— — / r{u) du-\- log r{B) 

u=0 ^ Jo 



1^2 hm 



— {u - u^)r{u) 



From (3.26) we obtain that r(0) = j3 and 

(/5 - f3')sn{uoo) 



lim (B - u^)r(B) = o ^ 
B-^uoo 2sn'(Moo) 



2cn(Moo)dn(tic 



hence 



1^2 



-2 lim 



B 



r{u) du — logr{B) 



-K, 



- (1 - OP 



-2-2 lim 



1 

— / r{u) du — log r{B) 
K Jo 



(3.54) 



(3.55) 



(3.56) 
(3.57) 

(3.58) 



(3.59) 



;i - OP. (3.60) 



(3.61) 



(3.62) 
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Using equation (3.50) for r{u), we immediately get that 



1 

K 







r[u)au — —— + log 



Ur. 



iK') 



(3.63) 



Now using equation (3.26) for r{u), we have 
lim 



1 

'k 



— ^+ hm log 



r{u) du — \ogr{B) 

Q{B + u^ + iK') (sn2(noo) - fin\B)) 



K 



K 



+ log 
+ log 



Q{B - Moo - iK'){fisn\u^) - p'sn^iB)) 



(3.64) 



= log 



Plugging this into (3.62) gives 



and thus we obtain that 



Z = -2 + 2 lo 



62 = 



2t?i(cu) 
261?! (a;)' 



(3.65) 
(3.66) 



4. RiEMANN-HlLBERT APPROACH: INTERPOLATION PROBLEM 



The Ricmann-Hilbert approach to discrete orthogonal polynomials is based on the follow- 
ing Interpolation Problem (IP), which was introduced in the paper [9] of Borodin and Bo- 
yarchenko under the name of the discrete Riemann-Hilbert problem. See also the monograph 
[5] of Baik, Kriecherbauer, McLaughlin, and Miller, in which it is called the Interpolation 
Problem. 

We will consider the lattice L„ defined in (1.27) and the weight Wn{x) defined in (1.25). 
Interpolation Problem. For a given n = 0, 1, . . ., find a 2 x 2 matrix- valued function 



l<i,i<2 



with the following properties: 



(1) Analyticity: Pn{z) is an analytic function of z for z e C \ L„. 

(2) Residues at poles: At each node x G the elements P„ii(z) and Pn2i(-2) of the 
matrix Pn{z) are analytic functions of z, and the elements Pni2(^) and Pn22(^) have 
a simple pole with the residues. 



Res Pnj2{z) = Wn{x)Pnjl{x), j = 1, 2. 
z=x 

(3) Asymptotics at infinity: There exists a function r{x) > on L„ such that 

0, 



lim r{x) 

x—*oo 



(4.1) 
(4.2) 
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and such that as z — > oo, Pn(-s) admits the asymptotic expansion, 



(4.3) 



where D{x, r{x)) denotes a disk of radius r{x) > centered at x and I is the identity 
matrix. 

It is not difficult to see (sec [9] and [5]) that the IP has a unique solution, which is 

p (~,\—( Pnn{z) C{WnPnn){z) \ / . . n 

where the Cauchy transformation C is defined by the formula, 

^(/)(-)=E7r^' (4-5) 

and Pnk{z) = + . . . are monic polynomials orthogonal with the weight Wn{x). Because of 
the orthogonality condition, as z ^ oo, 

oo / \ -p / \ oo oo j J oo 

x£Ln xGL„ j=Q j=k+2 

which justifies asymptotic expansion (4.3). We have that 

Kn = [Pl]l2, = [Pl]21- (4.7) 

5. Reduction of IP to RHP 

5.1. Preliminary Considerations. We would like to reduce the Interpolation Problem to 
a Riemann-Hilbert Problem (RHP). Introduce the function, 

27 . rnrz 

Il(z) = — sm — — . 5.1 



Observe that 



n(a;.) = 0, U\x,) = {-!)', exp (^^) = (-l)'^ for = ^ e 
Introduce the upper triangular matrices. 



(5.2) 



zmvz 



D:^(^)=(q ), (5.3) 

and the lower triangular matrices. 

Define the matrix-valued functions, 

rD"(^), when lmz>0, 
Rn = Pn(^)x|j^„^^^^ when lm;^<0, ^^"^^ 
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and 

f D^fz), when Im^; > 0, 
K-Pn{z)x{ -' (5.6) 

[ DL{z), when Imz <0. 

Prom (4.4) we have that 

[KUPn.-A^) - "■'''%r--''-' ^"'g +C-.C(.^..P„,„-i)(>-); • (5-7) 
when ± Im ^ > 0, 

and 

n(z)C(^/;„P„„)(^) 



n(^) Wn{z) 

"■'■Sr'"' - "■^-^'^fr-'W e-^'^ n(.)A-L.CKP.,„^OW ' (5-8) 



n(z) Uln(z) 

when ± Imz > 0. 



Observe that the functions Il'^{z), Rji(-2) are meromorphic on the closed quadrants of the 
complex plane and they are two-valued on the real and imaginary axes. Their possible poles 
are located on the lattice L„. An important result is that, in fact, due to some cancellations 
they do not have any poles at all. We have the following proposition. 

Proposition 5.1. The matrix-valued functions R"(-2) and R^iz) have no poles and on the 
real line they satisfy the following jump conditions at x 

(-1 nniwn{x) \ 
1 J' ^^-^^ 

and 

Rl4x)^RUx)fR{^), fR{^)-(_ii^ (5.10) 

Proof. It follows from the definition of R^(-2) that all possible poles of R^(-2) are located on 
the lattice L^- Let us show that the residue of all these poles is equal to zero. Consider any 
Xk e Ln. The residue of the matrix element R" ]^2(-^) is equal to 

Wn (•^fc)-fnn('^fc) / -j\k 

Similarly we get that 

Res RnM^) = 0, (5.12) 

Z=Xk 

hence R^i^) P^l^ x^- 

Similarly, the residue of the matrix element Rj^ ii{z) at Xk is equal to 

Kes K„,n(^J - -jZij-k ^^^^ - (5.13) 



Res RI,,{Z) = + Wr.{x,)Pnn{Xk) = 0. (5.11) 



In the same way we obtain that 



Res R„,2i(-2) = 0. (5.14) 

Z=Xk 
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In the entry Rj^ 21 (-^)) ^^e pole of the function C{wnPn)iz) si z — Xk is cancelled by the zero 
of the function Il{z), hence Rjj 2i(^) has no pole at Xk- Similarly, Rji 22('^) has no pole at Xk 
as well, hence RJj(-2) has no pole at Xk- 

Let us evaluate to jump matrices at a; G M. From (5.5) we have that 



which proves (5.9). Similarly, 

fj,{x) = T>i{xr'B^4x) 
which proves (5.10). 



■^2isin^ 

1 



mriwn{x) 



i 2i sin 1 

U{x)Wn{x) 27 





' 1 






mri 









(5.15) 

(5.16) 
□ 



5.2. Reduction of IP to RHP. Let us discuss how to reduce the Interpolation Problem 
to a Riemann-Hilbert Problem. We follow the work [5] with some modifications. Denote 



A = L„n[a',/3'], V = L„\A. 



(5.17) 



Consider the oriented contour E on the complex plane depicted on Fig. 7, in which the 
horizontal lines are Im^ = £,0, —e, where £ > is a small positive constant which will be 
determined later, and the vertical segments pass through the points z — a' and z — (3' . 
Consider the regions Jl^ and bounded by the contour E, see Fig. 7. Observe that the 
regions consist of two connected components, to the left and to the right of Q^. 



+ 



T 



a' 



1 



— ^ 


c 




— > 



Define 



Figure 7. The contour E. 

K^Rl(z)K-\ for zeQl, 
Rniz) = <; KMz)K-\ for z e nt 
KnP n{^)^n^, othcrwisc. 



where K„ 



1 

Q niri 

7 



(5.18) 



Define a contour T,r by adding to the contour E a vertical segment [is, —is], see Fig. 8. If 
A C C is a set on the complex plane and 6 e C then, as usual, we denote 



A + b ^ {z ^ a + b, ae A}. 



(5.19) 
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a' \ 


fO 



1 



- 









Figure 8. The contour T,r. 
Proposition 5.2. The matrix-valued function R„(^) has the following jumps on the contour 



R„+(2) = Rn-{z)jR{z), 



(5.20) 



where 



1 



when 2; e (— oo, a') U (/5', oo), 
(^_(l^)4;(^)-i 1) ^^^^ ze[a',P'], 



.0 



ij_ Wn(z)e 
nn n(2) 
1 



K,Di(^)K 



-1 



when z e (—00, a') U (/?', 00) ± is, 

n{z)-^ 

_|_ imrz 

mjL^^ n(z) 



(5.21) 



K„r>i{z)-'-Dl{z)K-' 



when z E (a , /3 ) ± i£ 
n(^) ;i^^n(^)e^ 2. 



ly/ien z G (0, iie) + a' or z E (0, iie) + [3\ 
^K„D°(;2)K;i w/ien ;2e(0,±i£), 



anc? 



D°(z) 



1 0' 



(5.22) 



Notice that the jumps on vertical contours close to the origin, D° (2;), are exponentially 
close to the identity matrix. 

6. First transformation of the RHP 
Define the matrix function T^n{z) as follows from the equation, 



Rn(^) 



-0"3 



T„(^)e 



"(9(2)-5)o'3 



(6.1) 
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where the Lagrange multipher I and the function g{z) are as described in Section 2 and 
(T3 = is the third Pauh matrix. Then Tn{z) satisfies the following Riemann-Hilbert 

Problem: 

(1) Tn{z) is analytic in C \ I^r. 

(2) T^+{z) = T^4z)jTiz) for z e Er, where 

gn(g_(^)-i)a3_^-^(^)g-n(5+(z)-l)a3 foj. 



(6.2) 



(3) 



Ti T2 

TJz)r^I-\ \ — - + ..., as z ^ 00. 

z z^ 



(6.3) 
(6.4) 



From (2.23) we have that 

g{z) —log z + 0{z~^), as 2; — > 00. 
This imphes that 

[Ti]i2 = e-"'[Ri]i2. (6.5) 
Let's take a closer look at the behavior of the jump matrix Jt described in (6.2) on the 
horizontal segments of Sr. We have that 

'^-nG(z) ^n{g+{z)+g-{z)-Viz)-ir 







:,nG{z) 



when z e (—00, a') U (/?', 00), 



Q-nG{z) Q \ 

.(^nny^-n{g+{z)+g-{z)-V{z)-l) ^nG{z) J whcu ze{a',(3'), 



1 ± 



„±nG(z) 



jriz) = < 



1-e 7 e T 



1 ± 



,n(2s(z)-i-V(z)) ' 



1-e T e T 



when z — X ±ie e (a, a') U (/?', /?) ± i£, 



when z — X ±ie & (—00, a) U (/3, 00) ± ie, 



(6.6) 



n(^)-i , 

H-e^^e-»(2sW-^(.)-^) n(^)) '^^^'^ ze{a,P)±ie. 



According to the properties of the g^-f unction, we have the following proposition: 
Proposition 6.1. The jump function jt has the following large n asymptotics: 

g-nG{z) g 
Q-nG{z) ]^ 



3t[z) = < 



gnG(z) 

1 0(e-"^(^) 

1 

1 e±'*^(^)0(e"~) 
1 



for ze{a',p'), 

for z e {a,a')U {(3',(3), 

for 2; e (— 00, a) U (/3, 00), 

for z e {a, a') U (/3', /3) ± ie, 



(6.7) 
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where C{z) is a positive continuous function on any subset of the given interval which is 
bounded away from the endpoints of each interval and satisfies 

C{z) > c\z + 1| for some c> 0. (6.8) 



7. Second transformation of the RHP 



The second transformation is based on two observations. The first is the well known 
"opening of the lenses" in a neighborhood of the unconstrained support of the equilibrium 
measure. Namely, notice that, for x e {a, a') U the jump matrix jri^) factorizes as 



-nG{z) 



1 

^nG{x) 

-- j-{x)jMj+{x 





-1 



1 

^-nG{x) ^ 



(7.1) 



which allows us to reduce the jump matrix Jt to the one Jm plus asymptotically small jumps 
on the lens boundaries. The second observation consists of two facts. Firstly, the jumps on 
the segments [a', (3'] ±ie behave, for large n, as ±e^^T~. Secondly, note that, for x E [a', (3'], 
G{x) is a linear function with slope With these facts in mind, we make the second 

transformation of the RHP. Let 

f T„(^)j+(^)-i for ^e{(Q;,a')U(/3',/3)}x(0,i£) 
Tn{z)j4z) for ze{{a,a')U{P',P)}x{0,-ie) 

imvz 

-^e 27 

nm 





Sn(^) 



7 



Tn(^) 

TJz) 



znivz 

^ 
otherwise. 
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for z e {a',p') X {0,ie) 
for z e {a',P') X {0,-ie) 



(7.2) 



This function satisfies a similar RHP to T, but jumps now occur on a new contour, E5, 
which is obtained from T^r by adding the two segments {a — ie,a + is) and (/? — ie, (3 + is), 
see Fig. 9. 









— ; 
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fa 

. — : 




: : 




p 



Figure 9. The contour E5. 
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On the horizontal segments for which the jump function js differs from Jt , we have that, 
as n — > oo. 

'or 

1 



_Q-nG(z) 



js{z) = < 



1 + 0{e~^^'^/'^) 0{e"'^~'-'^^^~^'^'^/'^^) 



^nG(z) 

1 + 0(e-^"^/^) 

rv7n^-n{2g{z)-l-V{z)) ^ _|_ (9(e~^"'^/'>') 

_g-n7ri(l+C) Q 
^-n{g+{z)+g-{z)-l-V{z)) _gn7ri(l+C) 



for z e {a,a)U{(3',p) 
for z-ie e {a, a') U /3) 
for z + i£ e (a, a') U /3) (7.3) 
for z e [a',f3'] ± ie 
for z 



By formula (2.34) for the G-function and the upper constraint (2.10) on the density p, we 
obtain that, for sufficiently small £ > 0, 



< TReGlx ± ie) = 27ip(x) + 0(s^) <— + 0(e'^' 

7 



(7.4) 



This, combined with property (2.32) of the (jf-function, imply that all jumps on horizontal 
segments are exponentially close to the identity matrix, provided that they are bounded 
away from the interval [a, For what follows we denote 



p3 

= TJ" + ^27r / p{x) dx — TT + n7r(l + (), 
Jo 



so that 



_g-n7ri(l+C) _ Q-i^n 



8. Model RHP 



(7.5) 
(7.6) 



The model RHP appears when we drop in the jump matrix js{z) the terms that vanish 
as n — > oo: 

(1) M{z) is analytic in C \ [a, ,8]. 

(2) M+{z) = M_{z)jM{z) for z e [a,/?], where 

r 

1 tor 2; G la.a'l U lU.U'l 

(8.1) 



jM{z) 



-1 

-io„a3 foj, ze[a',(3']. 



(3) 



, ^ Ml M2 

z z^ 



as 2; — > 00. 



(8.2) 



This model problem was first solved, in the general multi-cut case, in [12], and is solved 
in two steps. In the first step, we solve the following auxiliary RHP: 

(1) Q{z) is analytic in C \ [a, a'] U 

(2) Q4z) = Q_(z) for z E [a, a'] U [(3', (3] 
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(3) Q(z) = I + 0(z-^) as z ^ oo. 
This RHP has the unique solution (see [12]) 



/ 7(2)+7~M^) 7(^)-T~^(z) \ 
\ 2i 2 / 

where 

/ (.-«)(.-^-) y/^ (8 4) 

with cuts on [a, a'] U [/?', P]. 

To solve the model RHP described in (8.1) and (8.2), we again use elliptic functions. 
Define the function 

where is as defined in (1.38) with elliptic nome q = e*'^'^ = e 27 ( r = |^), and c? and c 
are arbitrary complex numbers. Notice that / has the following periodic properties: 

/(s + 7r) = /(s); f{s + nr)^e-'^J{s) (8.6) 

and that / is an even function. Now let 

where is as defined in (3.1). Then u is two- valued on and satisfies 

u^{x) — U-{x) = TIT for X G (8.8) 

Also, 

S±(«) = |, «±(aO = f±Y' ^±(/5') = ±Y' ^±(/^) = 0' (8-9) 
cf. Fig. 6. Because ^/R(x)^ — —y^R(x)_ for x E [a,a']U [P', (3], it immediately follows that 

u+{x) + u_{x) ^ for xe[/5',/5], (8.10) 

and that 

u+{x) + U-{x) = u+{a') — u+{j3') + U-{a') — U-{13') = t: for xE[a,a']. (S-H) 
We now define 



^3(^(^)+d+%^) ^3(-^(^)+c^+^k) 

1^3(^(2;) + d) ^?3(-'u(2;) + d) 

where d is an arbitrary complex number. It then follows from (8.6) and(8.8) that 

/i+(x) = e-^""/i_(x) and /2+(x) = 6^^-/2- (x) for xe[a',(3% 
and from (8.6), (8.10), and (8.11) that 

h+{x) = f2-{x) and f2+{x) = /i_(a;) for x e [a, a'] U 



8.12) 

8.13) 
8.14) 
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Define the matrix valued function 



1?3(u(2;)+d2) 1?3(-w(z)+(i2) 

where (ii and ^2 are yet undetermined complex constants. Then, from (8.13) and(8.14) we 
have that 



F^{x)^F4x)(^^~^^" for xe[a',P'] 

F+{x) = F_{x) (^^ ^ for X e [a, a'] U [/?',/?]. 

We can now combine (8.3) and (8.15) to obtain 



(8.16) 



iv±v^J iV^; I ^(^)_^-i(^)^?3(n(z)+d2 + %^) 7(^)+7-i(;;)i?3(-n(2)+d2 + %) I VO-^U 

2i i93(ii(z)+(i2) 2 i93(-w(z)+<i2) 



where 



■l93(Soc.+rfl + %) 



F(oo)=| ^3(..+..) .J. (8.18) 







'&3{-Uoo+d2) 

and Moo = u{oo). This matrix satisfies conditions (8.1) and (8.2) of the model RHP, but may 
not be analytic on C \ [q;,/3], as it may have some poles at the zeroes of 'd-s{±u{z) + di^2)- 
However, wc can choose the constants di and ^2 such that these zeroes coincide with the 

zeroes of 7(2;) ± ^^^(z) and are thus cancelled in the product. 

First consider the zeroes of '-/{z) ± '-/^^{z) These are the zeroes of 7^ (-2) ± 1 and thus of 
7^(2;) — 1, thus there is only one zero, which uniquely solves the equation 

which is 

-°= („,-^:;;g'-^.) "("'.'^'). (8.20) 

It is easy to check that 7(^0) = 1, thus xq is the unique zero of 7(2;) — 7~^(2;), whereas there 
are no zeroes of 7(2;) + ^~^{z) on the specified sheet. We use here the change of variables v 
defined in (3.2). Notice that, by (3.19), 

/ N P'-a dn^(Moo) .g^^v 



implying that 



2 



dn^f-u, 



sn u 



Since Xq G {a',f3'), we must have u{xq) G {iK', K + iK') (if we choose to take ■«+). Since 
sn^ is a one-to-one function on this interval there is a unique point uq G {iK\ K + iK') such 
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that sn^(xio) = p^^s^^^- The simple period identity 

sn(u + K + iK') = (8.23) 
kcYi[u) 

along with (8.22) gives that we must have 

Mo = u{xq) = K - Uoo + iK' (8.24) 

thus 

TT Tlx IX 

u{xQ) = —{K-u^ + iK') = — + --Uoo. (8.25) 
ZK Z Z 

We now consider zeroes of the function ^d^^u^z) — d) = ^d^^—u^z) + d). The zeroes of this 

function are the solutions to the equation 

TT TTT 

u{z)-d={2m + l)- + {2k+l)— (8.26) 

for any m,k E Z. Because u maps the first sheet of X to the rectangular domain [0, |] x 
[~"T' T"]' clear that this equation can have at most one solution, and without any loss 
of generality we may take m — k — 0. Then, if we want the solution of this equation to be 
xq, we need to let 

d ^ u{xo) - ^{1 + t) ^ -Uoo. (8.27) 

This choice of d also ensures that 'ds{u{z) + d) = '&3{—u{z) — d) has no zeroes on the first 
sheet of X. We can then let 

= d, d2^ -d (8.28) 

so that (8.17) and (8.18) become 



7(z)- 


2 

-7" 


^i(u(z)+d) 
Mu(z)-d+^) 


-2i 

7(z)+7-l 


1?3(-«(2)+rf) 
(^) ^3(-«{^)-d+%) 1 




2i 


'd-ci(u(z)-d) 


2 


1?3(-S(^)-d) / 






-1(2) ■ds.(u(z)+{n+\)Lo) 


7(z)-7- 




7(z)- 


2 
-7" 


'dA(u(z) + ^) 
j93(w(2)+(n-|V) 


-2) 
7(^)+7" 


; <94(«(z)-f) 

i?3(«W-(n-iM 




2i 




2 


^4{«(2)+|) y 



(8.29) 



M(^) = F(oo)-' 
= F(oo)-i 

where 

F(oo) = I ^^1"^ ^) I = I >^4(n.) I , (8.30) 

1?3(0) 

solving the model RHP. The asymptotics at infinity are 

M(^) = / + — + 0(^-2) (8.31) 
z 

where the matrix Mi has the form 





^93 (Woo +d+ )i?3 ( -«oo +d) 



* 



-4i 



-lyr _ I r'3{Uoc^-ta-t--^)rr3{-Uoo+a) i on\ 

^ ~ I ^3(^oc-ci+^)^3(-i»oc-ci) (/3-/3')+(a'-a) ^ I " ^^^""^^^ 

, 1?3 (-«oo -<i+ ^ )l?3 (Moo - Ci) 4i 

The matrix Mi can be written in a cleaner fashion and in terms of the original parameters 
as follows: 
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Proposition 8.1. We have that 



where 



For a proof of this proposition see Appendix D in the end of the paper. Notice that since 
M solves the model RHP, we have that 

detM(;s) = 1, zeC. (8.35) 



9. Parametrix at outer turning points 

We now consider small disks D{a,e) and D{f3,e) centered at the outer turning points. 
Denote D = D{a,e) U D{P,e). We will seek a local parametrix U„(^) defined on D such 

that 



(1) 
(2) 
(3) 



\Jn{z) is analytic on D \ E5. 
Vn+{z) = Vn-{^)js{z) for zeDHEg. 



V^{z) = M(^) (/ + 0{n-^)) uniformly for z e dD. 
We first construct the parametrix near (3. The jumps js are given by 

f 



js{z) = < 



^ ior z e {(3 - e,(3), 

for z e {/3,/3 + is), 



1 

_g-nG(z) ^ 



1 





for z e {13,13 — is), 



^-nG{z) ^n{g+{z)+g-(z)-V{z)-l) 
f,nG{z) 



for z e {(3,(3 \e). 



If we let 



then the jump conditions on Q„ become 

Qn+(^) = Q.n-{z)jQ{z) 



(9.1) 

(9.2) 
(9.3) 



(9.4) 



(9.5) 
(9.6) 
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where 



1 

-1 

1 

-1 1 

1 

1 1 
1 1 

1 



for z e {P- e,P), 
for z e {I3,l3 + ie), 
for z e {P,P — is), 
for z e (/?, /? + £). 



(9.7) 



where orientation is from left to right on horizontal contours, and down to up on vertical 
contours, according to Figure 7. 

Qn can be constructed using Airy functions. The Airy function solves the differential 
equation y" — zy, and has the following asymptotics at infinity: 



Ai(^) - 
M'{z) 



1 



_2 3/2 

e 3^ 



20F 



1/4 -2^3/2 

z ' e ^ 



1 + 



48' 



-3/2 



(9.8) 



as 2; — > oo with arg z e (— tt + £, tt — s) for any £ > 0. If we let 

yo(^) = Ai(2;), yx{z) = uK\{uz), y2{z) = uj^Ai{u;'^z) 
where u = e^, then the functions yo,yi, and y2 satisfy the relation 

yo{z) +yi{z) + y2{z) = 0. 



(9.9) 
(9.10) 



If we take 



yo[z) 



-y'2iz] 



for axgz e (0, 1) , 



-yi[z) -y2[z 

-y'liz) -y'^iz 

-y2{z) yi{z) 

-y'^iz) y[{z) 



for arg2; G ( ) ' 



TT 



for arg 2; e 



(9.11) 



then satisfies jump conditions similar to (9.7), but for jumps on rays emanating from 
the origin rather than from (3. We thus need to map the disk D{P,e) onto some convex 
neighborhood of the origin in order to take advantage of the function $^3. Our mapping 
should match the asymptotics of the Airy function in order to have the matching property 
(9.3). 

To this end notice that, by (2.19), for t e [/?', P], as t ^ 



pit) = C{P - tf'^ + 0((/5 - tfl^) , C > 0. 



(9.12) 
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It follows that, as z — > /3 for z e 



I 



P 2 
p{t)dt = Co(/3 - zf'^ + - zf'^) , Co = -C. (9.13) 



Thus 

V^M^) = -|y pW^tj (9.14) 

is analytic at /3, and so extends to a conformal map from D{(3,e) (for small enough e) onto 
a convex neighborhood of the origin. Furthermore, 

therefore ■(/'/g is real negative on [(3 — e,(3). and real positive on (/?,/5 + e). Also, we can 
slightly deform the vertical pieces of the contour E5 close to (3, so that 

e) n E5} = (-£, £) U (-i£, is). (9.16) 

Q„(z) = E^(^)$4n2/3VM^)) (9-17) 
U„(2) = E^(z)$^(n2/3V;^(2))e-"(^(^)-^-^)'^3 (9.18) 

E^(.) ^ Munur , Ku) - 5^ (""'""^f J (_\ :) . (9.19, 

where we take the branch of ip^J'^ which is positive on (/?, f3 + e) and has a cut on {/3 — 8,(3). 
We claim that E^(2;) is analj^ic in D{P,e), thus U{z) has the jump conditions of js- This 

is clear, as both M and have jump the same constant jump, (^^i ; the interval 

{(3 — s, (3], and are analytic elsewhere. The only other possible singularity for either M or 
is the isolated singularity at (3, and this is at most a fourth-root singularity, thus removable. 
It follows that 'Ei^{z) — 'M{z)Ij^{z)'~^ is analytic on D(P,e), thus U„ has the prescribed 
jumps in -D(/3, e). 

We are left only to prove the matching condition (9.3). Using (9.8), one can check that, 
for z in each of the sectors of analyticity, ^p{n'^^^ipp{z)) satisfies the following asymptotics 
as n ^ 00: 



We now set 
so that 
where 



1 A , M^y"^^ f-5 5t ^ ^, _2 



1 ^J^-^[-7 -7^)+^(-"^) 



X e 3' 

(9.20) 

where we always take the principal branch of ipj3{z)^/'^. As such, ipi3{z)^/'^ is two- valued for 
z e (/? — £, /3), so that 

= TTri / p{t)dt. (9.21) 

± Jx 
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Notice that, by (2.32) and (2.35), for x e {P - £,/3), 
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2g±{x) - V{x) = l±2m j p{t)dt. 

This imphes that 

[2^±(/?) - Vm - [2g±{x) - V{x)] = T27rz / p{t)dt . 

J X 

Combining these equations with (9.21) gives 



1 3/2 



± 



(25±(/3) - V{(5)) - {2g^{x) - V{x)) 



This equation can be extended into the upper and lower planes, respectively, giving 



{2g^{P) - V{P)) - {2g{z) - V{z)) 
Since, by (9.22), 2g±{(3) - V{f3) = I, we get that 



for ± Im z > 0. 



^^Jaizf'^ = -g{z) + ^ + 2 
for z throughout D{f3,e). Plugging (9.20) and (9.26) into (9.18), we get 

1 



(9.22) 



(9.23) 



(9.24) 



(9.25) 



(9.26) 



Vn{z) = M(z)L^(z)-i-^n-^-3^^(^)-^3 



1 A 

-1 i 48n 



-5 5i 
-7 -7i 



+ 0{n-^) 



^nig{z)-^-l)a3^-nig{z)-^-i)a3 



(9.27) 



= M(z) 



-3/2 



48n 



1 6i 
6i 



Thus we have that U„ satisfies conditions (9.1), (9.2), and (9.3). 

A similar construction gives the parametrix at the a. Namely, if we let 

2/3 



/ P{t)dt 

J a 



(9.28) 



then ipa is analytic throughout D{a, e), real valued on the real line, and has negative deriv- 
ative at a. Close to a, the jumps jg become 

1 1 



Q ^ J for 2; e (a — £, a), 

for z E {a,a + is), 
for z & {a,a — ie), 
for z & {a,a + e), 



1 
-1 1 

1 
1 1 



(9.29) 



1 
-1 
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where orientation is taken left to right on horizontal contours, and up to down on vertical 
contours according to Figure 7. After the change of variables ipa (and a slight deformation 
of vertical contours), these jumps become the following jumps close to the origin: 



1 

-1 

1 
1 1 

1 
-1 1 

1 1 
1 



for ipc^iz) e (-£,0), 
for i/jaiz) e (0,i£), 

for ipa{z) e (0, -ie), 
for iJa{z) e (0,£), 



(9.30) 



where orientation is taken right to left on horizontal contours, and down to up on vertical 
contours. These jump conditions are satisfied by the function 



1 
-1 



Then we can take 



U„(z) = E^(z)$4n^/=^V'a(^))e-"(^(^)-^-^)^^ 
for z e D{a,£), where 

E^(^) = Miz)K{z)-\ 



2v^ \ 



V4, 









n''',l^]I\z) 



-1 i 
1 i 



(9.31) 
(9.32) 

(9.33) 



Similar to (9.20), we have that in each sector of analyticity, ^a{'i^^^^'^a{z)) satisfies 



0-3 



1 -A i^ajz)-^" 

-1 -i) 48n 



-5 -5i 
-7 li 



+ Oin- 



(9.34) 

Once again, we have that, for x e (a, a + £), -0^(2^)^^^ takes limiting values from above and 
below, so that 

t)dt. (9.35) 





f 

= ±7ri 1 




± J a 



In analogue to (9.24), we have 

(2<,±(a) - Via)) - (2g(z) - V{z)) 



for ± Im 2; > 0. 



Since, by (9.22), 2g±{a) — V{a) = / ± vri, we get that 



\i^a{z) 



3/2 _ _ 



9{z) + 



V{z) I 



± m for ± Im 2; > 0. 



(9.36) 



(9.37) 
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Plugging (9.37) into (9.32) and (9.34) gives, as n — > oo, 
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+ 0(n- 



1 -A i'ajz) 

-1 -ij 48n 



U„(z) = M(^)L^(z)-i-^n5-3^^(^)-^3 



3/2 



-5 -5i 
-7 7i 



M(^) ^ 



-1 i 



-1 



.1 A^VO^fS 5z )^o(n-) 



1 i 



/ + 



V^q(2) 



-3/2 



48n 



-6i 



-6i 



48n 



+ 0(n- 



7 -7i 



(9.38) 



10. Parametrix at the inner turning points 

We now consider small disks D{a',£) and D{l3',e) centered at the inner turning points. 
Denote D — D{a',e) U D{P',e). We will seek a local parametrix lJn{z) defined on D such 
that 



(1) 

(2) 
(3) 



lJn{z) is analytic on D\T,s. 

Vn+{z) = Vn-{z)js{z) for zeDn^s- 



lJ„{z) = M(^) (/ + 0{n~^)) uniformly for z G dD. 
We first construct the parametrix near a'. Let 

U„(z) = Q„(z)e^^'^3e-"(^(")-^-5V3 for ±imz>0. 

Then the jumps for Q„ are 



1 
-1 

-1 

-1 -1 

1 -1 

1 

1 1 

1 



for z E {a' — e, a'), 

for z e (a', a' + e), 

for z e {a' , a' + ie), 

for z e {a', a' — ie), 



(10.1) 
(10.2) 
(10.3) 

(10.4) 



;io.5) 



where orientation is taken from left to right on horizontal contours, and down to up on 
vertical contours according to Figure 7. A proof of this statement is given in Appendix C. 
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yi{z) -y2{z) 



for axgz e 



(10.6) 



Then solves a RHP similar to that of Q„, but for jumps emanating from the origin 

rather than from a'. 

Notice that, by (2.19), for t e [a, a'], as i — > a'. 



p{t) = ^ - - ^)'^' + 0{{a' - tf'^), C > 0. 



It follows that, as z — > a', for z e (a, a'), 
1 
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- p(t) (it = C{a' - zf" + 0((q;' - zfl^) , Co = -C 



Thus, 



2/3 



(10.7) 



(10.8) 



(10.9) 



is analytic at a', and so extends to a conformal map from D{a\ e) onto a convex neighborhood 
of the origin. Furthermore, 

V'a'(a') = 0, C,(a')>0; (10.10) 

consequently, ■^q,/ is real negative on {a' — £, a'), and real positive on (a', a' + s). Again, we 
can slightly deform the vertical pieces of the contour E5 close to a', so that 

V'a'{^(Q;',£) nEs} = {-e,e)[J{-ie,ie) 



We now take 
where 



Q„(z)=E^'(z)$,,,(n^/=^Va'(^)) 
E^'(z) = M(z)e^^^3L„(z)-^ for ±lm^>0. 



20F \^ 







n 



1 i 
1 -i 



(10.11) 
(10.12) 

(10.13) 



and we take the branch of ip][i'^ which is positive on (a', a' + e) and has a cut on (a' — e, a'). 
Un then becomes 



\Jn{z) = M(^)e±'^'^^L^>)-i$«'(n'/'^a'(^))e^'^"^e-^(^(^)-'4^-5)-3 



2 2 - 
for ± \m.z > 0. 



(10.14) 
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The function {n'^^^'4'a'{z)) has the jumps js, and we claim that the prefactor E"' is analytic 
in D{a', e), thus does not change these jumps. This can be seen, as 



M+{z)e 



inn _ 



■0-3 



thus the jump for the function M(2;)e^ a"'^^ is 



e 2 "^jue 2 



e 2 



or equivalently, 



e^-3_^-^g^.3 



J^e"^"' for («'-£,«'), 

g^a3g-inna3g^a3 foj, ^ ^ ^^/^ + 



^ ^ ^ for 2; e (a' — £, a'), 

for 2; e (a', a' + £), 



(10.15) 



10.16) 



1 
1 



(10.17) 



which is exactly the same as the jump conditions for . Thus E" (2;) = M 
has no jumps in D{a', e). The only other possible singularity for E" is at a', and this singu- 
larity is at most a fourth root singularity, thus removable. Thus, E^ is analytic in D{a',e), 
and Qn has the prescribed jumps. 

We left check that U„ satisfies the matching condition (10.3). The large n asymptotics of 
$a/(n^/^'0a'(-2)) ^-re given in the different regions of analyticity as follows: 



+ 0(0 



i 1 



± 



T I 48n 
for ± Im 2; > 0, 



-5i 5 
7i 7 



;i0.18) 



where we always take the principal branch of ipa'izY^"^- As such, il>a'{zY^^ is two- valued for 
X & {a' — e, a), so that 



3/2 



± 



1 \ TVt 

p{t) j dt = =F— (a' - x) ±7ii I p{t)dt. 



Prom (2.32) and (2.35), we have that 

2g+{x) -V{x) = l + 2m p{t)dt , 2g_{x) - V{x) ^ I - 2m p{t)dt 

J X J X 

for X e (a' — £, a'). These equations imply that 

{2g±{x) - V{x)) - {2g±{a') - V{a')) = ±2m T p{t)dt. 

J X 

We can therefore write (10.19) as 



3/2 



± 



m , , s 1 
^2^(«-^) + 2 



{2g^{x) - V{x)) - {2g^{a') - V{a')) 



;io.i9) 



;io.2o) 



(10.21) 



(10.22) 
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We can extend these equations into the upper and lower half-plane, respectively, obtaining 



a - ^) + - 



{2g{z)-V{z))-{2g^{a')-V{a')) 



i / p{t)dt 

Ja' 



Using (10.20) at a; = a', we can write 

2 , . .0/2 Tri , , > , > V(z) I . i"^ 

or equivalently, 

2 , , , s Viz) I niz i(Qn — 7^) 

=»W-^-2±27=f-4;r^ for ±lm.>0. 

Plugging (10.18) and (10.24) into (10.14) gives 

Un{z) =M{z)e^'^'''L'^{z)-^^^n-'^''''i)a'{zy ' 



for ± Im z > . 

(10.23) 



for ±lmz>0, (10.24) 



(10.25) 



tT3 



X 



i 1 
i -1 



± 



' 2v^ 
48n 



-5i 5 
7i 7 



+ 0(n-2) 



X e 



M{z)e^'^'"'L^ {z)-^^n-'^'''ij^> (z)-^-3 

2vvr 



X 



1 t\ MzT^ 

1 -iy 48n 



5 -5i 
-7 -7i 



0{n 



-2\ 



= M(;2) 

= M(z) 



/ + 



-3/2 



48n 



_1 6.^^^.^.3 + 0(^-2) 



48n 



We can make a similar construction near Let 



2/3 



for ± Im (z) > 0. 

(10.26) 

(10.27) 



This function is analytic in D(/3', e) and has negative derivative at thus Im z and Im ipp'iz) 
have opposite signs for z e -D(/5', s). Then the jumps for Q„ are 



z = < 



1 
-1 

-1 

-1 -1 
1 -1 

1 

1 1 

1 



for ze{p',p' + e), 

for ze{p'-e,p'), 

for z e + is), 

for z e -ie), 



(10.28) 
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where the contour is oriented from left to right on horizontal segments and up to down on 
vertical segments according to Figure 7. After a slight deformation of the vertical contours 
and the change of variables ipfji, these jumps become the following jumps close to the origin: 



Jq(#'(^)) 



1 
-1 

-1 
-1 -1 

1 -1 

1 

1 1 

1 



for ip,3,{z) e {-e,0), 
for i>f}'{z) e {0,e), 
for ipi3'{z) e {—is, 0), 
for iljp'iz) e {0,ie), 



(10.29) 



where the contour is oriented from right to left on horizontal segments and down to up on 
vertical segments. These jump conditions are satisfied by the function 



;io.3o) 



Then we can take for z e D{P', e), 



for ± Im 2; > 0, 



(10.31) 



where 



J_fn-'/%y\z) 
20F 



n 



-1 
-1 



;i0.32) 



We once again have 



^f3'{n''/'iPf,.{z)) = ^n-"e-^iPp,{z)-"^''' 



I ~1 
^1. 1 



MzT^fbt 5 

48n \-7i 7 



+ 0(n- 



^■in^g,{zfr'a3 fQj, ^ Im^^/(^) > (so ± Im^ < 0), 



;i0.33) 



and for z e D{P',e), 



3V(^) = ±^ + ^(^)-^-2^^;^ for ±lm.>0. 



(10.34) 
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Combining (10.31), (10.33), and (10.34) gives 



X 



± 



20F 

1\ ^ ipp'iz)-^/^ (hi 5 



i 1 



48n 



7 + 0{n-') 



V(z) I . 



m{z)e 



/ + 



48n 



6i 1 



+ 0(n-2 

^ + 48n ive^e^'"^" 1 ) 



for ± \m.z > 0. 



11. The third and final transformation of the RHP 

We now consider the contour Ex, which consists of the circles dD{a, e), dD{a', e), dD{f3', e), 
and dD{P, e), all oriented counterclockwise, together with the parts of S5 \ ([a, a'] U [/?', /?]) 
which lie outside of the disks D{a,£), D{a' ,e), D{j3' ,e), and D{l3,e), see Fig. 10. 




Figure 10. The contour E^. 



We let 
X„(2;) = 



Sn(^)M(^)-^ for z outside the disks £>(«,£), D{a',e), D{p',s), D{(3,e), 

Sn{z)Vniz)-^ for z inside the disks £)(«,£), D{a',e), D{f3',e), D{f3. e). 

(11.1) 

Then X„(2;) solves the following RHP: 

(1) X„(2;) is analytic on C \ S^. 

(2) X„(2;) has the jump properties 

Xn+{x)=y.n-{^)jx{z) (11.2) 

where 

. 'M.(z)\Jn(z)~^ for z on the circles 

ix(^)- . ' . . (11-3) 



(3) As 2; — > 00, 



M{z)jsM{z)~^ otherwise. 



X„(^)~I+^ + ^ + ... (11.4) 
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Additionally, we have that jx{z) is uniformly close to the identity in the following sense: 

{/ + 0(n^^) uniformly on the circles 
(11-5) 
/ + 0(e-^(")") on the rest of Sx, 

where C{z) is a positive, continuous function satisfying (6.8). If we set 

fx{z)^jx{z)-I, (11.6) 

then (11.5) becomes 

{0(n^^) uniformly on the circles 
0(e"^(")") on the rest of Ex. 
The solution to the RHP for X„ is based on the following lemma: 
Lemma 11.1. Suppose v(z) is a function on Ex solving the equation 

.(.)=/- J- / "MM^du forzeEx (11.8) 
27^^ z_-u 

where z_ means the value of the integral on the minus side of Ex • Then 

X^{z)^I-^[ ^^^^^^dx. /or^eC\Ex (11.9) 
27rz Js^ z-u 

solves the RHP for X„. 

The proof of this lemma is immediate from the jump property of the Cauchy transform. 
By assumption 

X„„(^) =t;(z) (11.10) 

and the additive jump of the Cauchy transform gives 

X„+(z) - X„_(^) = v(z)f^(z) = X„_(z)ji(z), (11.11) 

thus X„_|_(2;) = X.n-{z)jxiz). Asymptotics at infinity are given by (11.9). 

The solution to equation (11.8) is given by a series of perturbation theory. Namely, the 
solution is 

oo 

v(z)^I + J2M^) (11-12) 

k=l 

where 

'^k{z)^-—. du, vo{z)^I. (11.13) 

'^'^'^ JT.X Z — U 

This function clearly solves (11-8) provided the series converges, which is does, for sufficiently 
large n. Indeed, by (11.5), 

\vkiz)\ < ( — j _^ , , for some constant C > (11-14) 
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thus the series (11.12) is dominated by a convergent geometric series and thus converges 
absolutely. This in turn gives 

oo 

X„(z) = 7 + J]X„,fe(^) (11.15) 

k=l 

where 

1 f Vk-l{u)j^{u) 



We will need to compute 



X„,i(^) = / ^^du. (11.17) 



12. Evaluation of Xi 

We are interested in the matrix Xi, which gives the ^-term of X„(2;) at infinity, see (11-4). 
By (11.9), 

Xi = -^ /" v(u)fAu)du, (12.1) 

hence by (11.12), (11.14), 



Xi 



27ri 

We would hke to evaluate the integral, 

1 



/ fx(u)du + 0{n-'). (12.2) 



27ri 





J 

s 



fx{u)du (12.3) 



with an error of the order of n^^. By (11.7), it is enough to evaluate this integral over the 
circles dD{a,e), dD{a',e), dD{P',e), and dD{(3,e). As we will see in the next section, the 
matrix- valued function j%{z) is analytic in the punctured disks, hence 



Xi = 



Res + Res 



Res 

Z=f}' 



Res 

z=l3 



3l{z) + 0{n-'). 



(12.4) 
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We will be especially interested in evaluation of the [12] element of the matrix Xi, and we 
will prove the following asymptotic formula. Introduce the numbers, 



^'^{nco + f ) 




VP' = - 



VP 



and 



MnuJ + l) Ml) + 




^4^ + 1)^3(1) 



2V2' 



-24 



^;(na; + fK(f) 



-24 



■i?3(m^ + 1 )^4(f ) 
(12.5) 



7.^, , 3,^ , 3, , , (a' - a)(/5 - a) 
-(/?'-«) + -(/?-«) + - (a' -a)- ^^,1'^^ ^ 

7,^ ,, 3,^, 3, , , (a' - a)(p' — a') 
-^{P -«)-2(/^ + ^^.-^^ , 



{(5 - a') 



Introduce also the numbers, 



(12.6) 



and 



^i(0)^i(nu; + f) 



Ml) ^4(na; + f) 

m) 



Ml) Mr>^ + i)' 



MO)^liru^ + l) 

m)^iir^) 

MO)^l{nu; + ^) 
Mi)^l{nu;) ' 



Ca' — 



m) 



Ml) MnoJ + ^) 
M^ + 1) 



m) 



^4(1) Mr^ + i) 



(12.7) 



(12.8) 



Lemma 12.1. ^45 n 



■ 00, 

[Xi]i2 



n 



iX^ + X^,+Xf,,+Xp) + 0{n-^), 



(12.9) 
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2{p' -a)J ' 



96 

96 



2{(3~ a') 



2(/5' - a) 7 ' 



(12.10) 



^ 96 V 2(/3-a') 
Proof of this lemma is given in the next section. 



13. Proof of Lemma 12.1 
On the circles, dD{a,e) and dD{P,e), we have 

^t^^^Miz)(^\ -_^I^M-\z) + 0{n-') for zedD{a,el 

r\z) + 0(n-^) for zedD(P,s) 



48n ^^^^6^ !;jM-^(z) + 0(0 for zeai5(/3,6), 

(13.1) 

and on the circles, dD{a',e) and dD{P',e), we have 

for ;2 e dD{a',e), ±lmz > 0, 
for z e dD(p',6), ±Imz > 0. 



Thus 



+ 0(0 

for z e dD{P',6), ±Imz > 0. 

(13.2) 



-"t^^^Miz) M-^(^) + 0(n-^) for . G aD(a,s), 



48n — -1 y ^ ^ lor ^ fc 

- ^'^48r^' M(^) M-^(;.) + 0(n-^) for z G aD(/5, 

- 48„ ^'"> (-WeT"- 1 )M--\z) + 0(,„-')CovzedD(a',sl 



48n ^ ^ 



/ -1 6ie±*^"\ , ^, _ 



for z e 9r>(/3', 



(13.3) 
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for ±Im 2; > 0. To simplify notation, we will write the model solution given in (8.29) and 
(8.30) as 



where 



^ ^3(^(^)-^^ + ^|.)^3(0) . ._ ^3(^(^)+ ^00-^)^3(0) 

''^'^ MHz) - uM^) MHz) + u^)M^) 

Muiz)+u^ + ^)Mo) MHz)-u^-^)Mo) 



(13.5) 



^3(m(^) + «oo)^3(^) ^ ' ^3(«(;2) - «oo)^?3(^) 



2 y '^jv"'v^; "'oo/^csv 2 

Notice that each of the functions -^ij is analytic throughout the complex plane, except on 
the intervals {a, a') and (/?',/3), where they satisfy the relations 

W=F, M± = (13.6) 
and on the interval (q;',/3'), where they satisfy 

[^n]+ - e-^^-[^n]- , [^^12]+ = e^''"[i?i2]- , [^?2i]+ = e-^^"[^?2i]- , [^22]+ = e^^'-M- 

(13.7) 

Multiplying out equations (13.3) gives 

Jxiz) = ^^^^^^f^(i' ^k)+0{n-') iorzedD{^,e), ^ = «,«',/?', /?, (13.8) 

^O'^ \J21 J22/ 

where 

Ju = I [3((f (^) + r\z)) i^i, - di,) + (f (^) - 7-^(^))^ii^i2 + - ^D] , 

JT2 = ^ [3(7^(z) +7-^(^))J(^) - (7^(^) -7-'(^))^ii^i2 + 6X(^)] , 



& = [3(7^(^) + r'{z)) J{z) + {j\z) - 7-'(^))^ii^i2 + 6K{z)] , 



(13.9) 



A = -2 [3((7'(^) + l-'{z)) {^i, - ^L) - h'iz) - r'{z))^nA2 + H^'n - ^D] , 
and 

f ^??ie^^" + e-^^" for Im^ > 0, 
J{z)^{ f .„ 'I .„ (13.10) 

X(^) = J (13.11) 
\ i??ie-*"" - i??2e' " for Imz < 0. 

In order to integrate jx{z), let us examine the behavior of the various functions described 
above near each of the turning points. Introduce the numbers 

Ac, = V(a' - «)(/?' - a){(3 - a) , A^, = ^/{a' - a){f3' - a'){(3 - a') , 

Ap, = ^{P'-a){P'-a'){P-P'), A^ = ^{(3-a){(3-a'){P-P'), ^^^'^^^ 



48 

and 
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B. 



1 1 1 

+ -77. h 



° a' — a (5' — a (5 — 

1 



1 1 

+ 



(3' -a (3' -a' (5 - (5' 
For a; e (/?,/? + £), we have 



Ba^ 



1 1 1 

+ -77. : + 



a' — a 13' — a' (3 — a' 
1 1 1 

+ -7^ 7 + 



(3 -a (3 -a' (3 - (3' 



u[x 



for X e (a — £, a), 



u[x) 



for X e {a', a' + e), 



u±ix) = ^u±{x) = I ± ^ - -^Vx - a' + 0{{x - a'f); 



and for x e {P' - e,P'), 



TTT TT 



2K 

Also, from (8.4) we have that 



An, 



{z - Ci){z - p') ^ l{z-a'){z-p) 





{a — 




3/2 


rj'\z) = 




a) 


-3/2 








-3/2 








-3/2 



(13.13) 



(13.14) 



(13.15) 



(13.16) 



(13.17) 



(13.18) 



{z-a'){z-P) \l {z - a){z - (3') ' 
and from (9.21), (9.28), (10.9), (10.27) , and (2.19) that 

^ + ^A^B^{a-z) + 0{{a-zr) 
^-^A^,B^iz-a') + 0{{z-af) 
^'^A,,B,,{P'-z)+0{{(3'-zr) 
^ + ^A,B,{z-P) + 0{{z-pr) 

It follows that the functions (7^ ± 'j~'^)^~^^'^(z) are meromorphic in a neighborhood of each 
of the turning points. In particular, at 2; = a, we have 



(13.19) 



(7^±7->.-^/^(^)=± 



{a' — a){P — a) 



2{(x-zf 



+ 



\(\ — z] 



{13' -a) ^ 3 

± — ( 3 ~ a) 

2 10^'^ ^ 



,3 I {a' - a){l3 - a) 

± — [a — a) ^ — : 

10^ ^^5 {(3' -a) 



(13.20) 



+ 0(1); 
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at z — (3, we have 



+ 



3_ _ I{l3-a){l3- 13') 



2{z-(3f (z-P) 



±^(/5-«') + ^(/5-«) 



10 



5 (/3 - a') 



(13.21) 



+ 0(1); 



at 2; = a', we have 

and at z — P', we have 

(f ±7-W/'(^) = ± 



{a' - a){P' - a') ^ 1 



2(z - a'Y (z - a') 



-10'°-°) + 5 (/a-) 



(13.22) 



+ 0(1); 



+ 



(/3' - ^) 



(/3^ -c.0(/3-/30 
2(/3' - zf 
^3 l (/3--a0(/3-/?0 
+ 10^^ ^^ 5 /3'-q; 



2 10^^ ^ 



+ 0(1). 



(13.23) 



Notice also, from the relations (13.5), that the functions -df^ +'&i2 and 'J?ii'i?i2 have no jumps 
in neighborhoods of a or f3, and take finite values at 2; = a and z = f3, thus are analytic 
in neighborhoods of a and (3. Using (13.14) and (13.15), we see that these functions have 
Taylor expansions about z = P, 



^l,{z) + ^Uz) ^ 2 



^i(O) ^i(Moo - ^) £ 



^?ii(;^)^?i2(;^) 



and about z — a, 



^i(O) r)Ua^ -J TT- a- / 
^i(I|^) ^i(^oo) ^2A|d^2l^^i 



3V"' 00; 



2^2 (^^2 



(^-/3) + 



2 j2 



rf2 / 



12 



(^-/?) + 



(13.24) 



^ii(^)+^?2(^) = 2 



+ 



^i(o) t^K^oo - ^: 



+ 



^11^12 



{a - z) + 
'a — z) + --- . 



(13.25) 

By a similar argument, J{z) and 'i?ii'i?i2 are also analytic in neighborhoods of a' and (3' 
and using (13.16) and (13.17) we can write their Taylor expansions about z — a', 



J{z) = 2 



00 " 2 ' 



rf2 



^i(%^) ^K^^oo) 



+ 



+ 



2^2, dM2 

7r2 rf2 



^2(0) ^2(^^ _ i 
^i(^) ^?(«oo) ' 2^2^rf«2 



^11^12 



[z-a')^ , 

{z - a') + 



(13.26) 
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and about z — P', 



J{z) = 2 



^l(O) ^i(Moo 7T^ d 



+ 



+ 



2Al, dv? 
2AI, dv? 



- ^) + 



(13.27) 



iP' -z) + 



Finally, notice that the function — -(9^2 is an odd function of u, and using (13.14) and 
(13.15), we can write, for x G + e), 



X 



47r ^^(0) ^l{u^ ( ^'^{u^) ^'^{u^ - ^) 



and for x G (o; — £, a), 

'47r ^liQ) ^1(moo-^) 



^3(1^00) i?3(iioo - %) 



+ 0(x - /3) 



(13.28) 



X 



^^4(^00) ^ 



Similarly, using (13.16) and (13.17), we can write, for x e (a', a' + s), 
K[x) = —\/x — a' 



+ 0(q; - x) 



(13.29) 



X 



and for x e {(3' - e,f3'), 
K{x) = yjl3' -x 



2 

2 



+ 0{x - a')) 



(13.30) 



X 



^2(«oo) ^'2(^00-^) 



+ 0(/3' - x) 



(13.31) 



Prom equations (13.19), (13.28), (13.29), (13.30), and (13.31), it follows that the functions 

(■^11 (^) ~ '^12 (^)) '^a^^^(^) and {i^li{z) — 'i?i2(^)) '^p^^'^i^) are meromorphic in neighborhoods 
of a and respectively, and have simple poles at 2; = a and z — (3, respectively, and that 
the functions K(z)iIj~^^^{z) and K{z)il)'^f^'^ [z) arc meromorphic in neighborhoods of a' and 
respectively, and have simple poles dX z = a' and z = j3\ respectively. 
Let us compute the residues of functions that appear in (13.9). Observe that 
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Prom (13.20), 13.21), (13.24), and (13.25), we obtain that 

_^i(0)^l(r^ + f; 



m' - a) - ^(,3 - a) - -(a' - a) + ^^'^^^ 



( 



and 



Z — p 



^i(0)^i(na; + f) 
^l(f)^l(na;) 



+ 



A{p-a') dv? 



Also, 



Res (^n(^)^i2(^)) h'{z) - r'iz)) ^-'/'{z) 



^i(Q)^l(r^ + f ) 



09' -a 3,^ 

■— H (3 - a) 

2 10^^ ' 



10 



{a' — a) — 



1 {a' — a){(5 — a) 
5 {(3' - a) 



+ 



4(/3' - a) 



and 



Res (^n(^)^i2(^)) {l\z) - r\z)) i^-^^^z) 

Z—p 



^l(0)^l(m^ + f) 



1 {(3 - a){(3 - 
5 (/3 - a') 



+ 



4(/3-q;') ( 



From (13.19), (13.28), and (13.29), we obtain that 

t?i(0)^9l(ncu + f ) 



Res6{^l,{z)-^l,{z))ij-'/\z) = -127r- 
and 

Res6(^?i(z) - = -127r 



^'M^ + f ) 



,^3(1) ^4(na; + 
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We now turn our attention to the residues at the inner turning points. Prom (13.22), (13.23), 
(13.26), and (13.27), we have 



Res^?,i,-J'\z)J{z){^\z)+r\z)) 

z=a' 



3(/?-a) + -(/3 -a)--{a - a) + j-^, 



(13.39) 



+ 



4(/5-a') dv? 



and 



Kes?,rp?'\z)J{z){^\z)+r'{z)) 

Z — p 



-3(/3 -a) + -(/5-/3)-- 



(13.40) 



37r2 rf2 



Also, 



4(/3'-q;) du^ 



z=l}' 



^^(0)^f(ncu + f) 
dl{^)§l{nuo) 



+ 



^2 



4(/3-q;') dv? 



{(3 - a') 



^11^12 



+ —{(5' - a') - —{a' - a) + - 



l{a' -a){l3' -a') 



10 



10~ 



13 -a' 



and 



(13.41) 



Res V^^(^)(f (z) - 7-^(^))^ii(^)^i2(^) 

Z — p 



-3/2, 
V \ 



+ 



d'' 



{(3' - a) 



^11^12 



+ {(3' -a')--{(3-(3') + - 



l{(3' ~a'){(3-(3') 



10 



10 



(3' -a 



z=l3' 



A{(3'-a) dv? 
Prom (13.19), (13.30), (13.31), we have 



and 



ResQ'il)~f''^{z)K{z) = 127r 



ResQ'il)pf'^{z)K{z) = 127r 



i^l(0)^i(no; + f) 



2\2. 



^^Uncu + I) 



^?2 



^9i(na; + |) 
t?i(f) i?2(na; + |) 



(13.42) 
(13.43) 

(13.44) 
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Combining (13.8), (13.9), (13.33), (13.35), and (13.37), we get that 



^?i(0)^l(na; + f) 



-127r 



^i(f)^l(m^) 

^3(1) Mr^ + ^) 



7,^, , 3,^^ , 3, , , (a' - a)(p - a) 



(3' -a 



TT 



4(/3' - a) dv? 



3(^?ll+^?12)- ^11^12 



Similarly, we have 



ResbK^)] = ^ 



•^i(0)^f(no; + f) 



^(/3-a') + ^(/3'-a')-^(a'-a) + 



(13.45) 
(«'-«)(/?'-«') 



+127r 



^2(^1 



t?i(ncu + |) 



TT 



Ci2 



4(/? - a') 



3J(z)-^?n^?i2 



Resbx(^) 

z=p' 



I 

96 



^l(0)^i(m. + f ) 



^(/?'-a) + ^(/?'-a')-^(/9-/?') + 



-127r 



and 



ResbK^)] 



96 



M%) ^2(na; + f) 



+ 



TT 



d2 



4(/3' - a) dv? 



(13.46) 
(3' -a 



(13.47) 



2>J{Z) - ^iii?i2 









4V 2 

/t^i(f) ^;,(na; + f) 



-^P-a)--{(5-a)--{(5-(5) + j--, 



+127r 

Using MAPLE for calculations, we get 



^2 



4(/? - a') dv? 



3(^n + ^?2)-^iA2 



^[3(^?i + ^y-^n^i2] 



= 2 



i?2(0)i?2(^eU + 



+ 7 



^U^^ + f ) 
^^inw + I) 



+ 17 



*^4(na; + f) ^3(1) 
^-(na; + |)^-3(f) 
i?4(na; + |)i?3(f). 



-24 



(13.48) 



(13.49) 



12 



^^9i(0)^9^(no^ + f) 



+ 7 



i^^(m^ + I) 
^^{nuj + I) 



+ 17 



4l 2' 



^4(1) 



^3(na; + |) <'(|) 
^3(m^ + f) ^4(1) 



-24 



(13.50) 



-^[3J(z)-^?n^i2] 



+ 7 



l?2(0)^?2(^^ ^ 1^ 



^2(|)^2(^^) 



.^^^(ncu + f) 



-5 



i?i(mi; + I) 



+ 17 



^2(1) 



^i(na; + f) ^21 2, 
2 



-24 



^^;(na; + |)^-(!) 
^?i(na; + |)i?2(f). 



(13.51) 
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32 



z=P' 



+ 17 



-24 



^'2(no; + f)t?;(f) 



(13.52) 



These formulae, combined with (13.13), prove Lemma 12.1. 

14. Large n asymptotic formula for hn 

We evaluate the large n asymptotic behavior of /i„„ and then we use formula (1.29). By 
(4.7), Kn = [Pi]i2, and by (5.18), 



[Pi] 12 = [R 



hence 



hnn — [R-l. 



1 12 



12 



nm 



7 



Furthermore, from (6.5) we obtain that 



12 



7 



nm 



7 



and from (7.2), that 

It follows from (11.1) that 
By (8.33), 



hnn — e"'[Si 



12 



mn 



7 



Si = Mi + Xi. 



[M 



lJl2 



'i?4(na;) 



vr(l + C) 



and by (12.9), 



where 



[Xi]i. = ^ + 0(n-), 



n 



c(n) = + + Xfj, + Xp 
is an explicit quasi-periodic function of n. Therefore, 



h = p 



nl 



iAd,{{n +l)u;) ^ c(n) ^ -2)' 



By (2.38), 



62 



7r^9;(0) _ A 
2e'diiuj) ~ 7 



7ri^;(0) 
2^?i(u;) ' 



7 



(14.1) 
(14.2) 

(14.3) 

(14.4) 
(14.5) 

(14.6) 

(14.7) 
(14.8) 

(14.9) 
(14.10) 
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hence 



2n 



'd4{nuj) n 

miA-"+'-ih{{n + 1)^') / ci(n) 
7e2""j?4(ria;) \ n 



nm 



7 



+ 0(n-2) 



where 



ci(n) 



c(n)'i94(nu;) 



iAi^4{{n+ l)uj) 
From (1.29) and the Stirhng formula we obtain that 

,2r^^ / _ \ 2n 



(n!)2(27)2'» V27/ 27rn 



1 



hence by (14.11), 



hr 



2n 



1 nvrA" 



(n!)2 V27/ 27rn 7e2»'(?4(na;) 

'^74(7^a;) \ n 



n 6n 



+ 0{n-') 



where 



G 



C2(n) = ci{n) - 



27 47i?i(w)' 
Observe that Ci(n) has the form, 

ci(n) = f{nw,u;), 

where f{x,uj) is a real analytic function, periodic with respect to both x and 
TT and 27r, respectively, so that 

f{x + TT, u;) = /(x, u;), /(x, u; + 27r) = /(x, a;). 

We can now summarize now the asymptotic formula for /i„/(n!)^. 

Proposition 14.1. As n ^ 00, 



n! 



_ Q2n+ 



,^4((n+l)a;) 



where 



and 



G 



1 + + 



n 



7r^;(0) 
471?! (a;) 



/o(a:,u;) = /(x,u;) 



6 



is a real analytic function which satisfies the periodicity conditions 

fo{x + 7r,uj) = fo{x,uj), fo{x,UJ + 2Tl) = fo{x,uj). 
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By (14.12), 

nn.,u)^^^^^l^^. (14.22) 

where 

X = X„ + X„/ + Xp, + (14.23) 
and explicit expressions for X^, Xa', Xpi, Xp are given in Lemma 12.1. 

In the subsequent sections, we carry out a concrete evaluation of f{nuj,uj), obtaining that 

/o(x,a;) = 0, (14.24) 

thus we can improve Proposition 14.1 to the following. 

Proposition 14.2. As n ^ oo, 

v% = G-'M^^±M(l+0(„-)), (14.26) 

where G is defined in (14-19). 

To this end, we will first show that f{nu!,u!) does not depend on n. 

15. f{x,Uj) IS CONSTANT IN X 

Denote ^ ^ 

z^nu^-, f{x,u;) = f(x - -,u;) , (15.1) 

so that 

f{z,u;)^f{nuj,Lo). (15.2) 
To prove that f{x,uj) is constant in x we will prove the following lemmas. 

Lemma 15.1. The function f{z,uj) is doubly periodic in z. 

Lemma 15.2. The function f{z,uj) is analytic throughout the z-plane. 

From these two lemmas, it follows immediately that f{z,uj) is constant in z, as it is a 
doubly periodic entire function of z, and it thus follows that f{x,u!) is constant in x. 
To prove Lemma 15.1, we will check that f{z + ttt, a;) = f{z,uj). By (14.12), (14.8), 

ci(n) = f{nw, u) = f{z, u;) = + + 1^3 + n, (15.3) 

where 

A"„i^i(//^') '^i(^')'ls(0)'^4('-) (n . lo^t . 



iA§^{nuj + u) 487ri?i(0)i?i(f )i?4(-z - ^)^^{z + |) V " 2(/3' -a))' 



ai?4(ria; + c^) 487ri?i(0)T?i(f )i?4(-z - f )i?4(-z + |) V " 2(/3 - a') / ' 

X^,t^4(na;) _ ^^(o;).^! (0)^1(2:) / , vr^r/^, 

' lA^^in^ + c^) 487ri?l(0)i?f (f )^?4(;^ - %)^^{z + |) V ^ ^' 2(/3' - a] 



iA-diinuj + u;) 487r'(?i(0)^?|(f )^?4(^ - + |) V 2(/? - a' 



;i5.4) 
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Observe that Yi, ^2, ^3, ^ can be written in the form 

Yi = Qiihu{z) + Qi2hi2{z) + Qishsiz) + Quhu{z), 

Y2 = Q2lh2l{z) + Q22h22{z) + Q23^23(^) + Q24^24(^), 
Y3 = Qsihsiiz) + Q32h32{z) + Qsahssiz) + Q34h34.{z), 

Ya = Q4ihAi{z) + QA2hi2{z) + QA3hAz{z) + Quhuiz), 

where 
hiiiz 



(15.5) 



hiA{z 

h23{z 
^32 

hAi{z 

h^A^Z 

and the numbers Qij do not depend on z. More specifically, 





/il2 




^'4(^)^4(^) 


hniz) 








^?4(^-f)?94(^+f)' 


dA{^-%)d,{z + ^y 




/i21 






h22{z) 


^[{z)Mz) 


^4(^-1)^4(^+1)' 






^4(^-1)^4(^ + 1)' 


^\{zf 


h2A 






h3i{z) 


^liz) 


^?4(^-f)^4(^+f)' 
















hsAiz) 




^4(-2-|)^4(-2+f)' 






^?4(^-f)^4(^+f)' 




hA2 






hAsiz) 




^4(^-1)^4(^+1)' 






M^-i)Mz+^y 


^'i{z)Uz) 












^4(^-1)^4(^+1)' 













(15.6) 



Q 



11 



487r^?;(0)^?i(|) 



a + 127r^^ + 



TT 



-5- 



^3(1) 2(/?'-a) I ^3(1 



+ 17 




Q12 
Qu 



4^1(0)^1(1) 4(/3' - a)^;(0)^i(f ) ' 
57n?i(a;)T?2(0) 



13 



96{P' - a)^UO)^i(l) ' 



96(/3'-a)^U0)^i(f)' 



and similar formulae hold for other Qij. In particular, notice that 

7r^i(a;)^i(0) ^1(^)^1(1) 



(15.7) 



Q13 + Ql4 
Q23 + <524 
Q33 + Q34 
Q43 + Q44 



8{P' - a)^; (0)^1(1) 8i)[{0)M^)M^)M^) 



8iP - a')^'i(0)^i(|) 8^U0)^i(f )^2(f )^3(|) ' 



(15.8) 



K/?'-a)^UO)^?(|) 
7r^?i(cj)^9i(0) 



^^'l(0)^l(f)^?2(f)i?4(|) ' 



8(/3 - aO^'i(o)^?l(|) s^[{o)M^)Mi)M^) ■ 
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Observe that all = Qij{oj) are periodic functions of uj of period 27r. It follows from 
equations (1.39) that the functions 

hj,{z), J = 1,2,3,4, (15.9) 

are doubly periodic, 

hji{z + tt) = hji{z) , hji{z + ttt) = hji{z), (15.10) 

while the functions, 

hj2{z), hj^{z), hji{z), J = 1,2,3,4 (15.11) 

satisfy the equations, 

h jk{z + 7r) = hjk{z), k = 2,3,4] 

hj2{z + -kt) = hj2{z) - 2ihji{z), hjs{z + ttt) = hjs{z) — Aihj2{z) — 4:hji{z), (15.12) 
h4{z + ttt) — hj4{z) — 4ihj2{z) — Ahji{z). 
This imphes that the functions Yj — Yj{z) for j = 1, 2, 3, 4, satisfy the equations 

Yj{z + n)^Yj(z), 

Yj{z + ttt) = Yj{z) + {-2iQj2 - ^Qjz - ^Qji)hji{z) (15.13) 
+ (-4.iQj3 - AiQj4)hj2{z). 
The proof of Lemma 15.1 then follows immediately from (15.13) and the following identities: 

(Qi3 + Qu)huiz) + (Q23 + Q2i)h2i{z) + (Q33 + Qzi)hzi{z) + (Q43 + QAi)h^i{z) = (15.14) 

(QlS + Qu)hl2{z) + (Q23 + Q24)h22iz) + {Q33 + Q34)h-i2iz) + (Q43 + Q44)h42{z) = (15.15) 
Quhuiz) + Q22h2l{z) + Q2.2hzi{z) + QA2hAl{z) = 0, 

which are proven below. Introduce here the notation 



The sum in (15.14) can be written as 

M^) [&l^l{z) - i^l^l{z) + ^l^l{z) - i^l^l{z)] 

which is zero by the Jacobi identity (E.9). 
The sum in (15.15) can be written as 



(15.16) 
(15.17) 



(15.18) 



(15.19) 



8^[{o)A^2^3^aM^ - ^)Mz + f ) 

Using the identities (E.2), we can write the expression in brackets in the numerator of (15.19) 



as 



^'Mz) - ^l^liz) - 4^l{z) + ^l^liz) 



+ 



Mz) 



(15.20) 



The first term in (15.20) vanishes by (E.9) and the second term vanishes by (E.IO). Thus 
(15.15) is proven. 



EXACT SOLUTION OF THE SIX- VERTEX MODEL 



59 



Finally, we can expand the sum in (15.16) and make the substitutions from identities 
(E.2), to obtain 



^1 



^3(0) + 



+ 



§1 



m) - 



mi 

1?2 



(15.21) 



Once again, the first term vanishes by (E.9). If we substitute the identity (E.ll) in the 
second and third terms, (15.21) becomes simply 

which is zero by (E.IO). This proves (15.16) and thus Lemma 15.1. 

We now turn to the proof of Lemma 15.2. Notice that in the fundamental rectangle 

(15.23) 



n^{z eC: ~ < Rez < ^ , Q<lmz < ttt]}, 



the function j{z^ lo) has the two simple poles, 

TTT 



^1,2 



2 2 



Because of the double-periodicity of /, we must have 

Res/(2;, a;) = — Res /(2;, a;). 



z=z\ 



Z=Z2 



Denote 

Then we have that 

§1 



Rij{oj) — Keshij{z,uj). 



R 



11 



, R 



12 



R 



21 



^'i(0)^i(^) 



, R22 — 



-^31 — Q/ /n\ Q / \ ' -^32 



R. 



■41 



^'i(0)i?i(a;) 



, R. 



■42 



(^1 - 2^4)^4 
^'1(0)^1 (a;) 

-(^^3-^^3)^3 



, Rl3 

, R23 



R 



14 



RoA — 



-"-33 — n/ /^\ n / \ ' -"-34 — 



(15.24) 

(15.25) 
(15.26) 

i&i - 2i'd'^ - ^4)^4 
^'i(0)^i(cc;) ' 



, R. 



■43 



, R. 



■44 



and 



-i&i - 2iA^ - ^2)^2 
(15.27) 

(15.28) 



Res f{z, a;) = QjkRjk- 

Z = Zl ' ' 

j,k=l 

A priori, the sum in (15.28) is quite complicated, so we evaluate first the imaginary part . 
Multiplying out (15.28) and again making the substitutions from (E.2), wc get 



Im Res/(^,a;) 

Z = Zl 



(15.29) 
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which is zero by (E.ll). 

Substituting the identities (E.5), along with (E.2), into the sum (15.28) gives 

„ ~, , 24:A + 17B + 12C + 10D + 7E + 5F 
Res J [z,uj) — 



Z=Zl 



(15.30) 



where 

B 
C 
D 
E 



^l^l^lmli^) - ^?^3(0) (^^^^(0) - ^1^)) 
-^3 (^^^1(0) +^^^^(0)), 

^^^3(0) {^l + ^4) + ^2^^ (^2(0) + <(0)) 



(15.31) 



mm) 



mM (^3 ^8 + ^i^^(o) {r, -^i+^t- 



Note that none of these terms involve derivatives of theta functions. We immediately have 
C = by the identity (E.8). We can also use this identity to write 



and (E.ll) to write 



mm¥m+^\m) 



(15.32) 



(15.33) 



(15.34) 



We now combine the terms A, 5, and E to obtain 

2AA + 175 + IE =2Ami^^l{^^im) - (^2^2(0) - ^4^4(0)) 

- 2i^^^^^^^^(o) + + ^1) (^^(0) + ^l(o)) 

-7^^(<^4(0)+^^^^(0)). 

By (E.ll) and (E.13) we can write this sum as 

2AA + 17S + 7£; = 2A^im^im){m2{^) - ^1^)) - 2M^i?^i?^i?^(0) 

+ 7^^^^(0) (^^ + ^t) - nl{m{^) + ^^^2(0)). 

We now combine all terms and use (E.ll) and (E.13) to write all terms solely in terms of 
'j?2,'»?4, and factors which are constant with respect to a;, yielding 



(15.35) 



2AA + 175 + lOD + 75 + 5F = 



m) m)m) 
x(m)+m)-m) 



+ 17 



^^^^(0) 

^i(o) 



(15.36) 



which is zero by (E.13). Lemma 15.2 is thus proven, and it follows that f{z,u}) is constant 
in z. To evaluate the constant, we can take z = 0. 
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16. Evaluation of /(O, o;) 

We will evaluate f{z,u) at 2; = 0. Notice that many of the functions hjk{z) vanish at 
= 0. In fact we have 



f{0,u) = Quhn{0)+Quhu{0)+Q23h23{0) + Q3^h3m+Q3,h3,{0) + Q,ih,^{0)+Quhu{0). 

(16.1) 

The identities (E.4) and (E.3) allow us eliminate all derivatives of theta functions from the 
sum (16.1) except 't?2>'*^2- Making these substitutions and simplifying, (16.1) becomes 



/(0,cu) 



+ ^1^3(0) - ^i^i(o)) / 5^'2'(o) i7(# 



where 



+ 



+ 



^i^h^h^'^l ^i^l'&3^l A^l'ds'&l 
'di{u){2AA + 175 + 7C + 5D + 3E + 2F) 



(16.2) 



c = m^imtmmM - ^1^2(0)), 
D = m [^^^^(0)^1(0) (^^^^(0) - + ^1(0) (^mmio) - «(o)) 

+mimim^M{o)+m3{Q)) 



(16.3) 



The identity (E.IO) imphes that all terms in (16.2) involving derivatives of theta functions in 
the numerator vanish. Additionally, (E.10)-(E.13) allow us to simplify the numbers B, C, D 
and E. Namely, we have 



B = 
C = 
E^ 



-mMmmi^i^ii^)+^i^im'-m3m)m)), 

-m.mmo), d = m2m2mi{o) (^m3{mi{o) + ^3^2(0)^^(0) + m2{omo) 

-m2m2mtm 

(16.4) 
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Combining these terms gives us 

2AA + 17B + 7C + 5D + 3E + 2F ^ 

Again using (E.10)-(E.13), this expression simphfies to 

24A + 175 + 7C + 5D + 3E + 2F = 8^l^l^l^l^l{0)^l{0)^l{0) . 
Inserting this into (16.2), we get 



/(0,W) = f{z,Uj) = f{7VjJ,Uj) 

by (E.l) and (E.6). It then follows that 



12i?;(0)i?ii?2i?3i?4 



1 



fo(nw,u;) = f{nw,u;) - g = 0. 
This proves (14.24) and therefore Proposition 14.2. 



17. Large n asymptotics of Z„ 
By substituting (14.25) into (1.24) we obtain that 



n-l 



2-n 



hi 



= 2»'ft„n[G-'M™(i + o(r^)) 



fc=l 



where C > does not depend on n. Thus, by (1.15), 
[sinh(7 — t) sinh(7 + i)]"'V, 



where 

F^2G sinh(7 - t) sinh(7 + t) 
Theorem 1.1 is proved. 



= Ci?4(na;)F" 1 + 0(n-^) , 



7rsinh(7 - t) sinh(7 + t)^'i{0) 
27^1 (c^) 
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Appendix A. Proof of formula (1.20) for t„ 



We have 



T„ = det 



l<i,j<n 



where 



sinh(27) 



2(1 - e^^l 



sinh(7 + t) sinh(7 - t) (1 - e-2(7-t))(l - e-2(7+0) 
2 2'^" 



l_e-2(7-t) 1 



-2(7+t) 



2 = 2^ e-2K7-t) + 2 ^ e-''(^+*) - 2 



i=0 



= 2^e-^'M + 25;^, 



-2«(7+t) 



i=0 



1=1 



Thus, 



■2j\l\ 
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(A.i; 



(A.2) 



(A.3) 



Prom equation (A.3) and multi-hnearity of the determinant function, we have that t„ is 
equal to 



2" J2 det 



^ g2Wi-27|h| ^2/i)e2*'i"2^l'il • 

(2/2)e2"2-27|«2| (^2/2)2e2*'2-27|/2| . 

(2Z3)2e2*^3-27K3| (2i3)3e2«3-27K3| . 



y(2/„)"-ie2*'"-27l'-l (2/„)V 



n„2ti„-27|«„| 



= 2" ^ A{21) Yl e2«i-2TKil ]^(2/,)^-^ 



i=l 



i=l 

n 



2"' 5^ J]A(/)e2«i-2^l'^in(^i)'"', 

;i,i2vi'n = — CO i=l 



(2Zi)"-ie2«i-27l'il ^^ 
(2Z2)"e2*''"^'^''^' 

(2/3)'^+^e2*^3-27K3| 



(•2/^)2n-2g2«„-27|Z„| 



(A.4) 



where A(/) = nt<j(^i ~ ^j) Vandermondc determinant. Note that, up to sign, this 

expression for r„ is invariant with respect to any permutation of I'^s. So, multiplying by their 
signs and then summing over all permutations, we get 



n 



= 2"^ E n e'"^"'""-' A(y Yl i-^r WiT^^'"^ 

h,l2,---,ln=—00 1=1 7re5n 



1=1 



thus 



J2 A(y'ne'"'"'^"''- 



(A.5) 



(A.6) 



ll,l2,---,ln =—00 
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Appendix B. Derivation of formula (1.24) 

Multilinearity of the determinant function, combined with the form of the Vandermonde 
matrix, allows us to replace A{1) with 



det 



/ 1 



1 

P2{l2) 



1 

Pl{h) 
P2{h) 



1 \ 

Pl{ln) 
P2{ln) 



(B.l) 



\Pn-l{h) Pn-lik) Pn-lik) ■■■ Pn-liQ/ 



where {Pj(a;)}°^Q is the system of monic polynomials orthogonal with respect to the weight 
w{l). Then (1.20) becomes 



2" 



ni 



2 OO 



(B.2) 



il,...,/n=l \7r65n 



k=l 



k=l 



The orthogonality condition ensures that, after summing, only diagonal terms are non-zero, 
so we get 



,2 OO 



n-1 



n! 



E E n p'ikyM n ^(^'^) = n 



(B.3) 



li,...,ln=l \7reS„ A;=l 



k=l 



k=0 



Appendix C. Proof of (10.5) 

From, (10.4), (7.1), (7.2), and (2.32), we have that the jump Jq on (a' — e, a') is given by 



g-n(g+(2)+9-(2)-V(^)-0\ 



^ '''''' ' .gn(g+(2)+5_(^)-y(^)-0 I e 27 

l' 
-1 



(C.l) 
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Prom, (10.4), (7.2), (6.2), and (5.21), we have that the jump Jq on (a', a' + e) is given by 



g— ^3g-n(g_(z) 5^-i)CT3^- gn(3+(z) 5^-5)<^3g — 



( inirz \ ^1 

— 27 \ / — -^e 27 

. in.. JT 
zmg-— y \^ -nEie— 

/ / \ V(z) I \ imrz _ 



^ / V 



riTri / \ 7 



rnvi 



^ ^n(g+{z)-^-l)a3 



v(,) fnni \ / 1 0\ / — ^ \ vi£) 



(C.2) 

Prom (10.4), (7.2), (6.2), (5.21), and analytic continuation of (2.32) into a neighborhood of 
[a, a'], we have that the jump jq on (a', a' + is) is given by 

^ g-"(9+W-^-5)<^3 /'-^ 0^ ^ gn(fl+W-|)<T3^-^g-n(5+(2)-f)a3^-_^^^-)-l 



/ / \ viz) I \ in-jvz _ 

X g"(9+(^) ^-2)'^3g^:^o-3 



gn^<T3 /"-^ 0^ ^ ^■^g-n(ff+(^)-i)a3^-^(^)-lgn(g+(z)-^-|)a3g^^3 







^ 












^ 






|(- 


— y 


1 



= g» — <^3 7 ^ ] I . '^T* 



7 7 



-n(g+(z)-i)a3 / ^ ^ ] en(s+{^)-^-|V3g^^3 



-2isin(^) -e-^^We'^A / e""^ 

gnn-)g^ g^ 1 I _gn(9+(^)+9-W-/-^) g"^ 



/„^~s ^(2:) V(^) iriTTz\ / v(z) \ 

-2isin(^)e" 2 _g-"^g— \ / e"" 2 \ --^^ 

-2ism(^) + e~e''^9+{z)+9-{z)-i-v{z)) \ ^ 

imrz imrz \ , 

g 27 _g 27 \ i-aixz /I —1 

e~ / VO 1 



(C.3) 
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Similarly, we have that the jump jg on {a' — ie, a') is given by 



imvz ^ ( I \ Viz) I I ( \ ^(2) l\ in-KZ _ 
j- — e^T^^e""*-^"^^'' 2^-2)'^3j^n{g-(z) 2-^-2)0-3g 2^0-3 




2i sin(^^)e"^ 

V(z) inivz 

2i sinf + e~~e"(^'+(^^+^-(^)~'~^(^)) 




Appendix D. Proof of Proposition 8.1 
Prom (3.16) and (8.7), we have 



TT 



(1-0 



This, combined with formula (7.5) for Vt^ immediately gives 

^z{uoo + d)^^{-u^ + d + %^) _ i?3(0)T?4((n + 1)0;) 



^z{-Uoo + d)§^{u^ + d + %^) 



Formulae (2.16) give that 

{(5 - (5') + {a' - a) ti^I^Q) 



Ai 4i 
Plugging the duplication formulae (E.6) and (E.7) into (D.3) yields 

iP - p') + {a' - a) nMO)MO)M<^) 



Ai 



2i 
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Combining (D.2) and (D.4), we can write the [12] entry of (8.32) as 

i7ri?4((n + IV) ^^3(0) MO)MO)M'^) 



mil, = 



2 i?4(na;) M^) M^) 



(D.5) 



2 Mnoo) ^?2($) Mn^) 
Similarly, we can write the [21] entry of (8.32) as 



W4((n — l)a;j 



Appendix E. Theta function identities 

The following identities (see [30]) are used in this paper. There are the identities involving 
derivatives of theta functions: 

^'l(0)=^2(0)i?3(0)i?4(0), (E.l) 



^4(^) 



Mz) - , (E.2) 



Mz) , (E.3) 



i»',{z)Mz) + Mz)Mz)^l(Q) 



- ^ 2.-.(^)M.)y.).l(o) ^ ^^.^^^^.^^^ ^ ^.^^^^.^^^^ 

= + 2.-2(.).3(j)y.).^(o) ^ MO) J^.^^^^.^^^ ^ ^.^^^^.^^^ j ^ 

^3(.) ^ »i . 2.,(.)..(;)y.).^(o) ^ ^ ^^.^^^^.^^^ ^ ^.^^^^.^^^^ 



(E.4) 
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the duplication formulae 



:e.5) 

^3(2;^)^?3(0)i?^(0) = ^i(;^)2^4(^)' + ^2(;^)'^3(;^)' (E.7) 
^4(2^)^1(0) = ^^(z) - ^t{z) = ^^(^) - ^t(^), (E.8) 

the addition formula 

My + ^)My - = ^l{y)^liz) + ^l{y)^l{z) 

and the identities relating squares of theta functions 



(E.9) 



i^l{z)i^l{0) = Mz)^liO) - ^2(^)^3(0), (E.IO) 
^^(^)^^(O) = ^^(^)^^(O) - ^l{z)^m, (E.ll) 
4iz)^m = M^mO) - ^l{z)m), (E.12) 

Mz)m) = - ^i{z)^m. (e.is) 
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